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ABSTRACT 


A  detailed  experimental  study  has  been  conducted  in  a  towing  tank  on  a 
3.048-meter-long  mathematical  model,  known  as  the  Wigley  hull,  to  study  the 
effects  of  viscosity  on  the  wavemaking  resistance  of  the  ship  form.  The 
measurements  included  total  resistance,  viscous  resistance,  pressure  distribu¬ 
tion,  and  boundary-layer  measurements  of  the  model  at  zero  trim  and  sinkage. 
The  three-dimensional  boundary-layer  measurements  extend  from  midship  to  a 
distance  of  0.1  model  lengths  downstream  of  the  stem  and  include  the  pressure 
distribution  on  the  body,  and  three  components  of  the  mean  velocity  measured 
by  means  of  a  five-hole  pitot  tube.  These  measurements  were  carried  out  for 
four  Froude  numbers  of  0.266,  0.313,  0.350  and  0.400. 

A  numerical  method,  using  the  small-crossflow  approximation  of  the 
boundary-layer  equations,  has  been  employed  to  calculate  the  boundary-layer 
characteristics  along  the  streamlines  on  the  hull  which  were  obtained  by  using 
Guilloton's  method.  A  comparison  between  calculation  and  experimental  results 
has  been  made. 

The  wave-resistance  coefficients  of  the  ship  form  have  been  calculated  by 
the  Michell  integral,  using  the  Michell  thin-ship  centerplane  distribution,  a 
hull  surface  distribution  and  a  slender-body  centerplane  distribution.  These 
calculations  have  been  carried  out  for  two  different  cases,  1)  wave  resistance 
calculation  without  the  effect  of  the  boundary  layer  and  wake,  and  2)  calcula¬ 
tions  including  these  effects.  These  calculations  have  been  carried  out  only 
to  first  order  and  have  been  compared  with  the  residuary  resistance  and  the 
wave-profile  resistance. 
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CHAPTER  I 
INTRODUCTION 

1 .  Background 

Since  Weinblum  (1930)  introduced  Michell's  theory  (1898)  to  the  field  of 
naval  architecture,  there  have  been  several  attempts  to  prove  its  feasibility 
by  comparing  the  computations  according  to  the  theory  with  measurements  in  the 
towing  tank.  If  one  plots  the  wave  resistance  of  a  thin  ship,  as  calculated 
by  Michell's  integral,  versus  Froude  number,  the  resulting  curve  has  pro¬ 
nounced  humps  and  hollows.  From  experiments,  however,  it  is  known  that  the 
measured  wave-resistance  curve  tends  to  be  much  smoother.  These  discrepancies 
between  the  measurements  and  theory  can  be  attributed  to  the  fact  that  the 
basic  assumptions  of  the  theory  are  not  fully  satisfied.  The  major  assump¬ 
tions  of  the  derivation  of  Michell's  integral  are  the  thinness  of  the  body, 
the  linearization  of  the  boundary  conditions  on  both  the  body  and  the  free 
surface,  and  the  inviscidness  of  the  fluid.  Since  water  is  viscous,  ship 
forms  are  usually  not  thin,  and  the  boundary  conditions  are  nonlinear,  these 
assumptions  are  approximate. 

In  the  last  decade,  with  the  development  of  computers  of  increased  speed 
and  capacity,  the  inviscid  problem  with  the  exact  hull  and  free-surface 
boundary  conditions  has  been  investigated  by  many  researchers.  Many  examples 
of  such  calculations  are  given  in  the  DTNSRDC  workshops  on  ship  wave- 
resistance  calculations  in  1979  and  1983,  and  the  International  Symposium  on 
wave  resistance  in  Japan,  1976;  and,  not  surprisingly,  the  results  are  in  much 
better  agreement  with  the  residuary  resistance  than  those  from  linearized 
theory .  This  has  led  some  researchers  to  suggest  that  wavemaking  is  indepen¬ 
dent  of  viscous  effects,  Newman  (1976),  and  that  viscous  effects  can  be  ne¬ 
glected  in  calculating  the  wavemaking  of  a  ship,  Gadd  (1976).  However,  these 
opinions  are  not  held  universally. 

There  are  strong  physical  reasons  for  believing  that  viscosity  plays  a 
role  in  wavemaking.  The  magnitude  of  crossflow  in  the  three-dimensional 
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boundary  layer  determines  whether  or  not  bilge  vortices  will  be  generated  near 
the  bow;  see  Tatinclaux  (1970).  Another  region,  where  vortex  formation  and 
separation  have  been  observed  by  Tzou  (1966)  and  Chow  (1967),  is  at  the  ap¬ 
proach  to  a  wave  crest  in  the  wave  profile  along  the  hull,  near  the  stern.  In 
both  cases,  the  pressure  field  about  the  hull  is  .modified  by  these  viscous 
phenomena,  so  that  the  wave  generated  must  also  be  affected.  In  addition,  the 
displacement  effect  of  the  boundary  layer  and  wake  modifies  the  outer  poten¬ 
tial  flow,  of  which  the  wave  potential  is  a  part.  This  effect  may  become 
especially  important  toward  the  stern  where  the  boundary  layer  is  thick,  and 
its  displacement  thickness  significant.  Havelock  (1948),  Wigley  (1963),  and 
Milgram  (1969)  have  shown  that  agreement  of  the  computed  wave  resistance  with 
the  residuary  resistance  can  be  definitely  improved  by  modifying  the  shape  of 
the  hull,  especially  near  the  stem,  in  accordance  with  the  displacement 
thickness  of  the  boundary  layer  and  wake,  and  applying  linearized  wave  theory 
to  the  so-thickened  body. 

Thus  it  has  been  found  that  agreement  with  the  experimentally-determined 
residuary  resistance  can  be  improved  either  by  calculating  the  wave  resistance 
more  exactly  by  inviscid-f low  theory,  or  by  introducing  viscous  corrections  to 
linearized  wave  theory.  In  an  attempt  to  shed  more  light  on  this  subject,  two 
investigations  were  undertaken  at  the  Iowa  Institute  of  Hydraulic  Research 
(IIHR).  In  the  first,  the  boundary  layer  and  wake  of  a  Series-60  ship  model 
were  thickened  by  roughening  the  hull  surface,  and  the  wave  resistance  was 
determined  by  means  of  surface  profile  measurements.  The  results  of  Moreno, 
et  al.  (1975),  and  Landweber  (1977)  showed  that  the  roughening  approximately 
doubled  the  viscous  resistance,  and  appreciably  decreased  the  wave  resistance 
by  40  percent  at  a  Froude  number  Fr  =  0.25,  by  50  percent  at  Fr  =  0.28,  by  13 
percent  at  Fr  =  0.31,  and  by  15  percent  at  Fr  =  0.34. 

In  a  second  investigation,  the  wave  resistance  of  a  thin  form,  for  which 
the  residuary  resistance  and  the  linearized  wave  resistance  had  been  calcu¬ 
lated  by  Weinblum  et  al.  (1952),  was  recalculated  in  various  ways  to  determine 
the  relative  importance  of  satisfying  more  exactly  the  boundary  conditions  on 
the  free  surface  and  the  hull  in  the  inviscid  theory,  and  of  including  the 
effect  of  the  boundary  layer  and  wake.  From  this  investigation,  Kang  (1978) 
concluded  that,  for  a  very  thin  form,  both  the  inviscid,  nonlinear  effects  and 
those  of  the  boundary  layer  and  wake  had  to  be  taken  into  account  to  obtain 
good  agreement  with  the  residuary  resistance. 
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In  the  present  study,  a  more  shiplike  form  has  been  selected.  The  se¬ 
lected  form  is  the  Wigley  parabolic  shape,  for  which  a  large  amount  of  exper¬ 
imental  data  and  computed  characteristics  are  available.  Shearer  (1965), 
Nagamatsu  (1979),  and  Adachi  and  Hinatsu  (1981).  There  are  other  advantages 
in  choosing  this  mathematical  form.  First,  the  boundary  conditions  on  the 
hull  surface  can  be  expressed  easily;  secondly,  the  Wigley  model  has  no  vortex 
system  shed  from  the  bottom  because  of  its  sharp  keel  and,  due  to  its  sharp 
bow,  there  is  no  wave-breaking  resistance;  see  Shahshahan  (1981).  Further¬ 
more,  it  has  been  shown,  Patel,  et  al.  (1983),  that  there  is  no  separation 
zone  at  the  stern.  Thus  the  total  resistance  may  be  considered  to  be  composed 
exclusively  of  wavemaking  resistance  and  viscous  resistance,  free  of  vortex 
formation  or  separation,  so  that  with  this  form,  attention  can  be  clearly 
focused  on  the  interaction  between  the  waves  around  the  hull  and  the  boundary 
layer. 


2.  Review  of  the  Literature 

Since  the  present  study  is  most  concerned  with  the  problem  of  how  ship 
waves  and  wave  resistance  are  affected  by  viscosity,  the  effect  of  the  nonlin¬ 
earity  of  the  boundary  conditions  will  not  be  discussed  here.  For  good 
reviews  of  nonlinearity  effects  refer  to  Newman  (1976),  and  Ogilive  (1976). 

Historically,  Havelock  (1935)  was  the  first  to  point  out  the  attenuation 
of  the  wave  due  to  viscosity.  He  defined  a  reduction  factor  to  the  wavemaking 
action  which  is  one  at  the  bow  of  the  ship  and  decreases  toward  the  stem. 
Another  method  for  viscosity  correction,  proposed  by  Wigley  (1937),  divided 
the  ship  wave  resistance  into  three  components,  due  to  the  bow  wave,  the  stern 
wave,  and  the  interference  between  the  bow  and  stern  waves.  He  defined  two 
correction  factors  for  viscosity  which  were  not  independent,  and  gave  empir¬ 
ical  expressions  which  defined  them  as  exponential  functions  of  Froude  num¬ 
ber.  Emerson  (1954)  considered,  instead  of  the  Wigley  assumption,  that  these 
two  factors  are  independent.  Furthermore,  he  employed  another  model  of  the 
viscosity  effect,  that  is  the  virtual  extension  of  the  stern.  A  similar  idea 
was  adopted  by  Inui  (1957)  but  with  a  slightly  different  form.  He  attributed 
one  of  the  reduction  factors  to  nonlinear  effects  of  the  finite  beam,  instead 
of  the  "thin  ship"  of  Michell's  theory.  He  called  it  the  self-interference 
coefficient  and  gave  a  semi -empirical  relation  for  it.  The  aforementioned 
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corrections  have  a  common  weaknesses,  i.e.,  they  lack  a  rational  theoretical 
background,  and  they  apply  only  to  certain  ship  forms,  and  lack  universal 
validity. 

The  earliest  analytical  treatment  of  the  problem  attributed  the  effect  of 
viscosity  to  the  displacement  of  the  flow  by  the  thickness  of  the  boundary 
layer.  It  was  done  for  the  first  time  in  the  work  by  Okabe  and  Jinnaka 
(1950),  and  a  practical  formula  was  given  by  Lavrentiev  (1951).  Since  the 
displacement  thickness  is  defined  only  at  the  body  surface,  the  addition  of 
this  displacement  to  the  hull  surface  results  in  a  body  which  does  not  close, 
since  the  displacement  thickness  continues  behind  the  body.  A  few  numerical 
examples,  given  by  Wigley  (1963),  showed  the  effect  to  be  too  small  to  account 
for  the  actual  phenomena.  He  then  considered  an  extension  of  the  boundary- 
layer  thickness  behind  the  body  in  such  a  way  that  the  virtual  streamline 
forms  a  closed  shape.  He  found  that  this  modification  had  a  remarkable  effect 
on  the  wave  resistance.  At  the  same  meeting,  Wu  (1963)  showed  another  numer¬ 
ical  example  but  the  result  did  not  agree  well  with  the  measured  values. 
Arguing  that  boundary-layer  effects  on  wave  formation  are  appreciable  over 
only  a  small  length  of  a  ship  form  near  the  stem,  Havelock  (1948)  suggested 
that  this  effect  be  represented  by  a  virtual  extension  of  the  hull  shape 
behind  the  stern  as  a  cusped  wake.  He  found  that  the  computed  wave  resistance 
with  this  modification  was  in  better  agreement  with  measurements.  Many  re¬ 
searchers,  following  Havelock's  lead,  suggested  some  sort  of  wake  model  behind 
the  body.  The  proposed  wake  models  are  mostly  based  on  the  inviscid  flow,  and 
may  be  classified  into  two  groups.  These  are  called  the  potential-wake  models 
and  rotational-wake  models  according  to  Maruo  (1976),  who  wrote  a  thorough 
review  of  these  methods.  Potential-wake  models  are  those  in  which  the  effect 
of  displacement  of  the  wake  is  represented  by  some  distribution  of  virtual 
sources  or  dipoles.  The  sources  that  represent  the  displacement  by  the  wake 
are  called  Betz  sources.  They  were  considered  first  by  Betz  (1925)  in  his 
formulation  for  the  experimental  determination  of  viscous  drag.  The  aforemen¬ 
tioned  Havelock  modification  of  the  stem  shape  corresponds  to  the  addition  of 
these  sources.  Since  the  tail  point  of  the  modified  shape  becomes  a  cusp, 
this  type  of  wake  is  called  a  cusped  wake.  It  is  unnecessary,  however,  that 
the  modified  source  distribution  generate  a  closed  form,  and  the  wake  may  be 
of  finite  width  up  to  infinity  downstream.  Milgram  (1969)  employed  such  an 
open  wake  model  to  calculate  the  wave  resistance  of  a  ship  form. 
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Rotational  wake  models  assume  some  distribution  of  vortices  in  the  wake 
region,  and  a  general  theory  for  the  wave  formation  by  such  a  vortex  distribu¬ 
tion  was  given  by  Brard  (1970).  Tatinclaux  (1970)  carried  out  a  computation 
of  the  wave  resistance  of  an  ogival  strut  accompanied  by  a  rotational  wake  by 
assuming  a  volume  distribution  of  vertical  line  vortices.  Beck  (1971),  as¬ 
sumed  another  type  of  rotational  wake  composed  of  plane  vortex  sheets . 

The  effects  of  viscosity  have  also  been  investigated  by  Lurye  (1968), 
Wilson  (1971),  Maruo  (1976),  and  Mori  (1978),  who  linearized  the  Navier-Stokes 
equations  of  the  wake  and  obtained  solutions  in  terms  of  Oseenlets  represent¬ 
ing  the  wake  flow.  However,  since  it  has  been  shown  by  Gadd  (1963),  and  Swain 
and  Landweber  (1974)  that  the  wavemaking  of  vorticity  is  negligible,  and  by 
Preston  (1954),  Lighthill  (1958),  and  Landweber  (1978,  1981)  that  the  irrota- 
tional  flow  equivalent  to  a  boundary  layer  and  wake  can  be  generated  by  either 
certain  distribution  of  sources  or  doublets,  or  as  the  flow  about  a  form 
thickened  by  the  displacement  thickness,  it  appears  that  the  more  sophisti¬ 
cated  and  difficult  calculations  employing  Oseenlets  would  not  yield  results 
for  wavemaking  essentially  different  from  those  employing  the  displacement 
thickness.  Here  "equivalent"  denotes  that  the  real  flow  and  the  "equivalent" 
irrotational  flow  coincide  exterior  to  the  boundary  layer  and  wake.  In  spite 
of  these  many  efforts,  a  satisfactory  solution  of  viscous  effects  on  ship 
wavemaking  is  not  yet  available. 

3.  The  Objectives  and  Method  of  Approach 

The  problem  to  be  studied  is  that  of  the  effect  of  the  boundary  layer  and 
wake  of  a  ship  on  its  wavemaking  resistance.  The  objectives  of  this  research 
are,  first,  to  determine  the  importance  of  viscous  effects  on  the  prediction 
of  ship  wavemaking  resistance,  and  second,  to  incorporate  viscous  effects  into 
the  procedures  for  computing  wavemaking  resistance  in  a  way  which  is  physi¬ 
cally  reasonable  and  compatible  with  available  procedures  for  inviscid  flow. 

In  the  present  study,  in  accordance  with  1979  David  W.  Taylor  Naval  Ship 
Research  &  Development  Center  ( DWRNSRDC )  workshop  recommendations,  the  Wigley 
ship  model  will  be  restrained  so  that  there  is  no  sinkage  or  trim,  and  its 
resistance  and  flow  characteristics  will  be  measured  in  a  towing  tank.  The 
measurements  include  total  resistance,  viscous  resistance  by  means  of  a  wake 
survey,  the  pressure  distribution  on  the  hull,  and  velocity  profiles  in  the 
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three-dimensional  turbulent  boundary  layer.  Using  these  data,  the  boundary- 
layer  characteristics  of  the  hull  will  be  calculated.  These  results  will  be 
compared  with  a  numerical  solution  of  the  three-dimensional  boundary  layer 
along  the  streamlines  on  the  hull.  Guilloton's  method  (1951)  will  be  employed 
to  obtain  the  streamlines  on  the  hull. 

The  next  step  will  be  the  calculation  of  the  wave  resistance  of  the 
model,  by  using  both  Michell's  thin-ship  and  a  slender  body  centerplane  source 
distribution.  In  this  part,  we  shall  consider  two  different  cases,  1)  wave- 
resistance  calculation  without  the  effect  of  the  boundary  layer  and  wake,  and 
2)  calculations  including  the  effect  of  the  boundary  layer  and  wake.  These 
calculations  will  be  carried  out  only  to  first  order  and  will  be  compared  with 
the  residuary  resistance. 


4.  Plan  of  the  Report 

Chapter  II  contains  a  description  of  experimental  details,  including  the 
instrumentation  and  data-acquisition  systems.  The  results  are  presented  and 
discussed  in  Chapter  III.  Numerical  procedures  for  boundary-layer  calcula¬ 
tions  and  a  comparison  of  the  results  of  measurements  with  the  calculated 
values  are  discussed  in  Chapter  IV.  In  Chapter  V,  waveraaking-resistance 
calculations  with  various  refinements  are  presented  and  compared  with  experi¬ 
mental  values.  Chapter  VI,  contains  the  conclusions  and  the  recommendations 
for  future  work. 
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CHAPTER  II 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


1.  The  Model 

All  experiments  were  performed  in  the  Iowa  Institute  of  Pfydraulic 
Research  (IIHR)  towing  tank  which  has  been  described  by  Martin  (1958).  The 
towing  tank  is  91.44-m  long,  3.048-m  wide  and  3.14-m  deep.  The  ship  model 
employed  in  this  study  was  a  member  of  the  series  of  mathematical  models  with 
parabolic  sections  and  waterlines,  introduced  by  Wigley  (1942).  This  form  is 
characterized  by  sharp  edges  at  the  bow,  stern  and  keel,  and  is  defined  by 


M  =  f  [1  -  (f21)2 ]  (1  -  (|)2] 


(II 


-  vf  <.  x  <_  ^  -  H_<zj<0 


in  a  Cartesian  (x,y,z)  coordinate  system  with  the  origin  at  midship  on  the 

undisturbed  free  surface,  increasing  value  of  x  opposite  to  the  direction  of 

the  ship's  motion,  and  the  z-axis  directed  vertically  upward.  Here  L  is  the 

length  of  model,  B  is  its  beam  and  H  the  draft.  The  selected  form  had  a 

length  of  3-048-m,  and  parametric  values  of  P  =  0.100  and  P  =  0.0625.  The 

I*  L 

block  coefficient  of  the  model  is  Cg  =  0*444  &nd  its  wetted  surface  area  is 
"-•381  m  .  The  model  was  equipped  with  three  rows  of  pressure  taps,  from 
midsection  to  the  stem  on  the  portside.  These  rows  were  at  the  levels  |z|  = 
0.635,  3.81,  and  11.43  cm.  For  turbulence  stimulation  along  the  hull,  a  row 
of  plastic  studs  of  3.2  mm  diameter,  1.6  mm  height  and  9.5  mm  spacing  was 
fitted  on  the  model  at  15.2  cm,  five  percent  of  the  model  length,  from  the 
bow.  Figure  1  shows  the  Wigley  model.  With  the  towing  arrangement  used,  the 
model  was  restrained  in  both  trim  and  sinkage  for  all  the  experiments,  except 
the  total-resistance  measurements,  which  were  performed  for  both  fixed  and 
free  conditions. 
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Figure  i.  Wigley  hull  model. 


2.  Instrumentation 


A  five-hole  yaw  probe  was  used  to  obtain  the  direction  and  magnitude  of 
the  velocity  in  the  three-dimensional  boundary  layer  and  wake.  The  probe  was 
L-shaped  with  a  shorter  portion  of  75  mm  and  a  longer  section  of  300  mm.  It 
was  made  of  five  stainless-steel  tubes  of  one  millimeter  inner  diameter,  with 
overall  diameter  of  4*2  mm.  Since  it  has  been  shown  by  Bryer  et  al.  (1955) 
that,  for  boundary-layer  measurements  in  general,  a  sharp-edged  probe  is  more 
satisfactory  than  a  rounded  one,  as  it  is  less  affected  by  changes  of  Reynolds 
number,  the  five— hole  probe  was  fabricated  with  a  sharp  tip.  For  moving  the 
probe,  a  traverse  mechanism,  which  could  move  the  probe  in  three  directions  of 
a  Cartesian  coordinate  system,  was  designed  and  fabricated  in  the  II HR  shop. 
The  lateral  and  vertical  movements  were  measured  by  counters  with  the  minimum 
counts  of  2.00  and  2.05  mm,  respectively.  The  pressure  tubes  from  the  five- 
hole  yaw  probe  were  connected  to  five  differential-type  pressure  transducers, 
which  were  mounted  on  a  rigid  fixture  on  the  traverse.  In  order  to  obtain  the 
optimum  sensitivity  in  measurements,  these  transducers  were  of  different 
ranges.  The  center  hole  of  the  probe  was  connected  to  a  pressure  transducer 
with  a  range  of  0.02  N/m2  (0.3  psid),  and  the  other  holes  were  connected  to 
pressure  transducers  with  a  range  of  0.01  N/m2  (0.15  psid).  These  pressure 
transducers  formed  a  component  of  a  strain-gauge  bridge.  The  out-of-balance 
voltages  from  these  bridges,  which  are  proportional  to  the  pressure  sensed  by 
transducers,  then  were  sampled,  digitized  and  recorded  by  the  data  acquisition 
system. 

In  the  total-resistance  measurements,  the  model  was  attached  to  the 
dynamometer  of  the  towing  carriage.  A  force  transducer  with  a  range  of  ±  4.54 
kgf  was  employed  to  measure  the  total  resistance  of  the  model  via  a  frequency 
counter. 


3«  Calibration 

3.1.  Force  Transducer 

The  force  transducer  was  calibrated  against  dead  weights.  This  calibra¬ 
tion  was  checked  before  and  after  each  set  of  total  resistance  measurements, 
in  orders  of  both  decreasing  and  increasing  weights.  The  calibration  was 
found  to  be  linear  and  very  stable. 
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3.2.  Pressure  Transducer 

The  pressure  'transducers  were  calibrated  against  an  alcohol  manometer 
with  a  least  count  of  0.025  mm  of  alcohol.  The  calibrations  were  performed  at 
the  experiment  site  using  the  actual  connecting  cables  and  data-acquisition 
system.  During  the  course  of  measurements,  three  or  four  points  of  the  cali¬ 
bration  curves  were  checked  every  other  day.  The  calibration  curves  were 
found  to  be  linear  and  very  stable. 


3.3.  Calibration  of  5-Hole  Yaw  Probe 


The  five-hole  yaw  Probe  was  calibrated  in  an  open- throat  wind  tunnel. 
The  probe  was  placed  in  the  working  section  of  the  tunnel  on  sin  assembly  that 
allowed  it  to  be  yawed  as  well  as  pitched  about  the  flow  direction  at  any 
prescribed  angle.  A  total-head  tube,  mounted  at  the  same  longitudinal  loca¬ 
tion  but  a  small  spanwise  distance  away  from  the  yaw  probe,  served  as  the 
reference.  The  five-hole  probe  was  connected  to  a  single  differential  pres¬ 
sure  transducer  through  a  scanivalve;  see  Cramer  (1984).  The  static  pressure 
tap  of  the  transducer  was  connected  to  the  static  pressure  of  the  total-head 
probe.  Since  it  has  been  shown  by  Rajaratnam  and  Muralidhar  (1968)  that  the 
effect  of  Reynolds  number  on  calibration  is  insignificant  at  air  speeds 
greater  than  12  m/sec.,  the  calibration  was  carried  out  at  an  air  speed  of 
17.554  meters  per  second. 

If  the  uniform  flow  is  assumed  to  be  along  the  x-axis,  the  yaw  angle  \p 
and  pitch  angle  0  are  defined  as  shown  in  Figure  2.  This  figure  also  shows 
the  numbering  scheme  used  for  the  five-hole  probe.  Let  P-j_,  ?2>  P3>  P4»  and  P5 
be  the  pressures  sensed  by  the  holes  numbered  one  through  five,  respec¬ 
tively.  Then  it  is  evident  that  the  holes  numbered  2  and  4  are  primarily 

sensitive  to  yaw  and  the  holes  numbered  3  and  5  are  sensitive  to  pitch.  Based 
on  this,  and  extending  the  concept  used  by  Ramaprian  et  al.  (1978)  with  a 
three-hole  yaw  probe,  the  following  coefficients  are  defined: 


K 


similarly  for  pitch  angle: 


K 


(P2+P4) 


+  V 


(2) 

(3) 
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Furthermore,  if  Q  is  the  magnitude  of  the  stream  velocity  and  Pg  the  local 
static  pressure  at  the  probe,  the  velocity  coefficient  is  defined  as 


P1  -  4  (P2+P3+P4+P5) 


2  p  Q 


(4) 


and  finally  the  static-presure  coefficient  is  given  by 


K 


S 


P  -P 
1  S 

1  2 

7  p  Q 


where,  in  the  calibration  process,  Pg  =  0. 


(5) 


The  probe  was  calibrated  for  the  range  of  -35°  _<  35°  and  -35°  <  0  < 
35°  at  5-degree  steps  in  yaw  and  pitch  angles.  The  aforementioned  nondimen- 
sional  calibration  coefficients  were  calculated  from  these  measurements.  The 
nondimens ional  yaw  coefficient  was  plotted  against  the  yaw  angle  >P  for 
different  pitch  angles.  Similarly,  K0  was  plotted  against  9  for  different  yaw 
angles.  Also  Kq  and  Kg  were  plotted  against  0  for  different  values  of  ip .  See 
figures  3  through  6. 


Later,  the  five-hole  yaw  probe  was  calibrated  dynamically  in  the  towing 
tank  at  zero  angles  of  yaw  and  pitch;  see  Figure  7. 


4.  Data  Acquisition  System 

At  the  heart  of  the  data-acquisition  system  is  an  HP-1000  E-Series  mini¬ 
computer  which  was  connected  to  the  experiment  site  through  a  16-channel 
Preston  analog-to-digital  converter.  The  carriage  velocity  and  pressure  data 
were  sampled  simultaneously  by  this  system.  A  block  diagram  of  the  system  is 
shown  in  Figure  8.  A  generator  linearly  converts  carriage  speed  to  frequency, 
which  is  then  transformed  into  voltage  by  a  frequency -voltage  converter.  The 
sensed  pressures  from  the  5-hole  yaw  probe  are  converted  to  voltages  by  pres¬ 
sure  transducers  and  amplified  by  signal  conditioners.  These  voltages  are 
sampled  simultaneously  through  six  different  channels  of  the  analog-to-digital 
converter,  digitized,  and  then  read  and  stored  by  the  computer. 
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Figure  2.  Definition  sketch  for  the  yaw  probe  (pitch  A  yaw) 


Figure  3.  i aw  probe  calibration  curves  for  pitch  angles. 
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Figure  4#  Yaw  probe  calibration  curves  for  yaw  angles. 
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Figure  5.  Yaw  probe  calibration  curves  for  the  magnitude  of  velocity.. 
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Dynamic  calibration  curves  for  five-hole  yaw  probe. 
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Figure  8.  Data-acquisition  system  for  the  carriage  speed  and 


5.  Experimental  Techniques 

5.1.  General  Remarks 

The  model  was  attached  under  the  carriage  to  the  dynamometer.  For  the 
free  condition,  there  was  only  one  point  of  attachment  and  two  guides  were 
used  at  the  how  and  stern  to  prevent  lateral  movements  of  the  model.  The 
model  alignments  were  checked  by  levels  and  also  by  measuring  the  distances  of 
the  bow  and  stern  from  the  wall  of  the  tank  to  make  sure  that  the  model  was 
not  yawed  or  pitched  with  respect  to  the  direction  of  motion.  For  the  fixed 
condition,  a  parallelogram  arrangement  was  used  instead  of  guides  at  the  bow 
and  stern  to  prevent  the  model  from  deviating  from  its  initial  sinkage  and 
trim. 


Since  a  system  of  water  and  air  was  used  to  transmit  pressure  to  the 
pressure  transducers,  there  was  the  problem  of  air-bubble  trapping  in  the 
connecting  tubes.  This  was  avoided  by  cleaning  the  tubes  before  each  experi¬ 
ment  by  compressed  air  and  then  adding  a  drop  of  photoflo  fluid  to  each  tube. 

The  Analog-to-Digital  (A/D)  converter  of  the  data-acquisition  system  can 
read  +  10  volts  and  it  is  recommended  that  the  entire  range  be  used  to  reduce 
cable  noise  at  low  voltages  caused  by  the  long  cable  length  of  210  meters 
connecting  the  experiment  site  to  the  A/D  converter.  For  this  reason,  a  gain 
of  600  was  used  to  magnify  the  low  voltages  from  the  pressure  transducers. 
Cable  noise  was  found  to  affect  the  measured  pressure  by  ±  0.57  x  10-5  kg/cm2 
in  the  first  four  channels,  and  by  ±  0.88  x  10“^  kg/cm2  in  the  next  two  chan¬ 
nels  of  the  A/D  converter.  The  relative  magnitude  of  these  noise  levels  to 
the  measured  quantities  was  found  to  be  less  than  0.1  percent. 

For  tests  in  the  towing  tank,  a  time  interval  of  15  minutes  between 
carriage  runs  was  found  to  be  necessary.  This  was  determined  by  injecting  dye 
beneath  the  free  surface  and  studying  the  motion  of  dye  filaments.  It  was 
observed  that  even  after  15  minutes  the  water  has  a  velocity  about  0.2 
cm/sec.  This  gives  an  error  of  0.12  percent  in  the  measured  velocity. 

5.2.  Total-Resistance  Measurements 

Initially,  the  mechanical  system  used  to  tow  the  model  introduced  large 
fluctuations  in  the  reading  of  data.  This  was  due  to  speed  oscillations  of 
the  speed  control  at  low  speeds,  which  resulted  in  inertial  forces  on  the 
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model  of  the  same  order  of  magnitude  as  the  resistance  to  be  measured.  .Some 
preliminary  tests  showed  that  these  oscillations  are  of  high  frequency,  so 
that,  by  using  a  low-pass  filter,  this  problem  was  corrected.  Recorded  read¬ 
ings  with  and  without  the  electronic  filter  showed  a  marked  improvement  of  the 
system.  For  resistance  measurements,  the  output  signal  of  the  transducer  was 
averaged  over  a  time  interval  of  10  or  5  seconds  depending  on  the  speed  of 
that  particular  run.  At  least  two  time  averages  of  the  signal  were  recorded 
during  each  run. 

Total  resistance  measurements  were  performed  for  a  range  of  Froude  num¬ 
bers  up  to  0.400,  for  free  and  fixed  conditions.  These  results  are  shown  in 
Tables  1  and  2. 

5.3.  Boundary  Layer  and  Pressure  Measurements 

Initially,  it  was  attempted  to  measure  the  pressure  distribution  and 
boundary  layer  on  the  Wigley  hull  via  a  scanivalve.  It  was  found,  however, 
that  the  duration  of  a  run  in  the  towing  tank  was  not  long  enough  for  pressure 
readings  of  a  five-hole  pitot  tube  to  attain  their  asymptotic  values.  Hence, 
the  scanivalve  was  abandoned,  and  five  pressure  transducers  were  directly 
connected  to  the  5-hole  pitot  tube.  Pressure  distribution  on  the  hull  was 
also  measured  by  direct  connection  of  transducers  to  the  pressure  taps.  With 
this  arrangement,  mean-flow  measurements  were  taken  in  ten  streamwise  loca¬ 
tions  from  midship  to  the  stern.  At  each  section  the  boundary-layer  measure¬ 
ments  were  carried  out  for  five  waterlines.  Figure  9  shows  the  measuring 
positions  on  the  Wigley  hull.  These  measurements  were  continued  into  the  wake 
of  the  model  for  four  streamwise  sections  and  five  girthwise  locations  at  each 
section.  The  boundary-layer  and  wake  measurements  were  made  at  four  different 
speeds,  corresponding  to  the  Froude  numbers  0.266,  0.313,  0.350,  and  0.400. 
For  each  Froude  number,  since  the  model  velocity  was  not  exactly  the  same  for 
different  runs,  the  pressure  readings  of  the  5-hole  pitot  was  corrected  to  the 
speed  corresponding  to  the  Froude  number  by  the  following  relation: 

V  2 

P  =  (^t)  P'  (m  =  1,5)  (6) 

m  V 1  m  7 

c 

where  Vf  is  the  carriage  speed  corresponding  to  the  measured  value  of  P!  and 
c  im 

v0  is  the  speed  corresponding  to  the  desired  Froude  number. 
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TOTAL  RESISTANCE  COEFFICIENT  Ct 
(FIXED  CONDITION  T-18.3  C) 


xaaaaaaaaaaaaasaasaaass 


FR 

Ct 

FR 

Ct 

0.100 

0.00440 

0.300 

0.00496 

0.120 

0.00434 

0.315 

0.00511 

0 .140 

0 . 00431 

0.330 

0.00503 

0.160 

0.00424 

0.340 

0 . 00500 

0.180 

0.00425 

0.350 

0.00500 

0  .  200 

0.00429 

0.360 

0.00500 

0.220 

0.00432 

0.370 

0.00505 

0.240 

0.00450 

0.380 

0.00520 

0.265 

0.00451 

0.390 

0.00540 

0.280 

0 . 00464 

0.400 

0.00560 

Table  1 . 

Total  resistance  coefficients  of 
(fixed  condition) . 

Wigley  model 
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TOTAL  RESISTANCE  COEFFICIENT  Ct 
(FREE  CONDITION  T=13.3  C  ) 


FR 

Ct 

FR 

Ct 

0.106 

0.00455 

0.285 

0.00500 

0.123 

0.00448 

0.303 

0.00538 

0.139 

0.00456 

0.314 

0.00553 

0.158 

0.00457 

0.331 

0.00538 

0.183 

0.00458 

0.344 

0.00528 

0.198 

0.00455 

0.360 

0.00537 

0.213 

0 .00447 

0.372 

0.00557 

0.233 

0.00457 

0.385 

0.00590 

0.250 

0.00473 

0.400 

0.00625 

0.273 

0.00477 

Table  2 


Total  resistance  coefficients  of  Wigley  model 
(free  condition). 


Ul 


Figure  9.  Body  plan  and  measuring  positions  of  Wigley  model 
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CHAPTER  III 

RESULTS  AND  DISCUSSION 


1.  Error  Analysis  of  Measurements 

Before  discussing  the  results  obtained  from  the  measurements  of  velocity 
field  and  pressure  distribution  around  the  hull,  it  may  be  necessary  to 
analyze  the  error  involved  in  the  measurements  in  order  to  have  a  better 
understanding  of  the  results.  Generally,  the  error  involved  in  the  exper¬ 
iments  can  be  divided  into  two  groups:  physical  errors,  which  are  encountered 
during  the  course  of  measurements,  and  numerical  errors,  which  occur  in  pro¬ 
cessing  the  resulting  data.  In  the  following,  a  brief  description  of  the 
sources  of  physical  and  numerical  errors,  together  with  their  order  of  magni¬ 
tude,  is  given. 

1.1.  Physical  Errors 

The  most  common  sources  of  error  in  measuring  either  dynamic  or  static 
pressure  in  the  present  study  are  the  following: 

i)  Alignment  error:  The  yaw  angle  of  the  5-hole  yaw  probe  was  checked 

against  a  level,  and  its  pitch  angle  was  checked  against  a  parallel 
line  by  eye.  This  may  give  about  a  one-degree  error  in  pitch-angle 
alignment.  Meriam  and  Spaulding  (1935),  have  shown  that  the  resulting 
error  in  measuring  dynamic  pressure,  for  one-degree  misalignment,  is 
very  much  less  than  one  percent  of  actual  value. 

ii)  The  influence  of  probe  support:  The  pressure  gradient  associated  with 
the  curvature  of  flow  lines  around  probe  supports  was  avoided  by  plac¬ 
ing  the  probe  at  least  about  15  diameters  from  its  support,  which  is 
three  times  the  recommended  five-diameter  distance.  Thus  this  source 
of  error  was  completely  avoided. 
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iii)  Wall-proximity  effects:  When  a  total-head  probe  approaches  a  wall,  it 
is  shown  by  MacMillen  (1956)  that  the  streamlines  are  deflected  by  the 
probe-wall  interaction.  Such  an  interaction  results  in  an  error  of  up 
to  10  percent  in  dynamic  pressure  measurement  when  the  distance  from 
the  axis  of  the  probe  to  the  wall  is  less  than  1.5  probe  diameters, 
i.e.,  in  the  present  measurements,  the  first  two  measured  points  of 
boundary  layer  at  each  location  are  affected  by  this  error. 


iv)  Pitot-tube  displacement  effect:  It  is  shown  experimentally  by  Young 

and  Mass  (1936),  and  theoretically  by  Hall  (1956),  that  a  pitot  tube 
whose  axis  is  at  a  fixed  point  A,  measures  the  total  pressure  at  a 
point  B  a  short  distance  away  from  A.  This  distance  which  is  known  as 
displacement  effect  of  the  pitot  tube  due  to  its  finite  size,  in  gen¬ 
eral  is  a  function  of  the  shape  of  the  tube,  its  diameter,  kinematic 
viscosity  and  velocity  of  the  fluid.  For  the  5-hole  probe  used  in  the 
present  study,  the  ratio  of  the  displacement  effect  to  the  probe  diam¬ 
eter  is  approximately  between  0.1  for  Fr  =  0.400  up  to  0.14  for  Fr  = 
0.266.  This  means  that  the  shift  in  measured  points  due  to  displace¬ 
ment  effect  of  five-hole  pitot  tube  is  less  than  0.5  mm. 

v)  Turbulence  effects:  Turbulence  may  affect  the  reading  of  the  pitot 

tube  in  two  ways.  Firstly,  the  r.m.s.  value  of  the  fluctuating  veloc¬ 
ity  components  may  directly  affect  the  calibration  of  the  pitot  tube. 
However,  this  is  generally  believed  to  be  small  and  of  secondary  impor¬ 
tance.  Of  primary  importance  is  the  direct  effect  of  turbulent 
velocity  components  on  the  total-head  tubes,  which  causes  the  sensed 
value  of  mean  total  pressure  to  increase.  This  can  be  shown  by  a 
simple  analysis  of  the  measurements  of  a  pitot  tube: 


P 


m 


1  —2 
P  +  j  p  Q 

n  1  2  1- 

p+2pq  +  2’pq 


(7) 


Unfortunately,  in  the  present  study,  the  turbulence  quantities  were  not 
measured  and  an  exact  estimate  of  the  turbulent  effect  is  not  pos¬ 
sible.  However,  it  has  been  shown  by  Goldstein  (1936),  and  Walsche  & 


27 


Garner  (I960),  that  the  pitot  tube  responds  mainly  to  the  longitudinal 
component  of  turbulent  velocity,  especially  when  the  resultant  mean 
velocity  does  not  make  too  large  a  pitch  or  yaw  angle  with  the  axis  of 

_  2 

probe.  So  one  can  simply  replace  the  total  turbulent  energy  term  q'  , 
by  the  longitudinal  component  u '  .  For  the  usual  turbulence  intensity 
encountered  in  the  water  flows  the  turbulence  effect  on  the  probe  is 
small;  e.g.,  for  a  turbulence  intensity  of  20  percent  which  is  too  high 
for  practical  situations,  the  pressure  reading  of  the  yaw  probe  will 
increase  by  only  2  percent. 

vi)  Hole  geometry  for  static-pressure  taps:  According  to  Shaw  (I960),  an 
'ideal'  tap  geometry  is  a  small  circular  hole  of  less  1/4-mm  diameter 
drilled  perpendicular  to  the  surface  on  which  the  pressure  is  to  be 
measured,  with  the  corner  of  the  hole  perfectly  sharp  and  squared 
off.  Any  departure  from  this  geometry  will  introduce  errors.  Errors 
with  a  practical-sized  hole  occur  because  of  flow  in  and  around  the 
hole  opening.  In  the  present  study,  the  diameter  of  the  pressure  taps 
was  one  millimeter.  This  should  introduce  an  error  of  less  than  one 
percent  of  -|pU^  compared  with  a  hole  of  mm.  Nonperpendicularity  of 
the  pressure  taps  with  the  wall  will  also  introduce  some  error,  but 
this  error  for  small  deviation  is  of  secondary  importance  and  can  be 
neglected.  Thus  the  overall-error  estimation  due  to  the  geometry  of 
the  pressure  taps  in  this  case  is  less  than  one  percent. 

1.2.  Numerical  Errors 

As  already  mentioned  in  Chapter  II,  a  computer  program  was  developed  in 
order  to  obtain  the  velocity  vector  in  the  boundary  layer  from  the  five-hole 
yaw-probe  measurements.  In  the  numerical  scheme  used,  polynomials  of  the 
third  degree  were  fitted  by  the  method  of  least  squares  to  the  calibration 
results,  and  values  at  desired  points  were  obtained  by  linear  interpolation. 
To  obtain  an  estimate  of  the  errors  due  to  the  polynomial  fittings  and  inter¬ 
polations,  a  numerical  study  was  conducted.  A  series  of  calibration  data  with 
zero  pressure  coefficients  and  exactly  known  pitch  and  yaw  angles  was  used  as 
an  input  to  the  program.  The  output  for  this  test  showed  that,  for  small 
angles  of  yaw  and  pitch  between  zero  and  ±  5  degrees,  the  error  was  less  than 
0.4  degree.  For  angles  between  ±  5  and  ±  15  degrees,  the  error  varied  between 
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0.5  to  1.2  degrees.  The  error  increased  to  2  degrees  for  pitch  or  yaw  angles 
between  15  and  25  degrees.  For  larger  pitch  and  yaw  angles  between  25  and  35 
degrees,  the  error  was  between  2.2  and  3*5  degrees.  These  numerical  errors 
are  indeed  very  small,  first  because,  for  most  part,  the  flow  around  the 
Wigley  hull  can  be  classified  as  one  with  small  cross-flow,  and  the  measure¬ 
ments  showed  that  the  yaw  or  pitch  angles  are  less  than  10  degrees  over  most 
of  the  hull,  except  near  the  stem  where  they  can  be  as  high  as  15  degrees. 
Secondly,  since  these  angles  are  used  in  trigonometric  functions  to  calculate 
the  velocity  components,  these  numerical  errors  become  even  less  appreci¬ 
able.  This  test  showed  that  the  overall  numerical  error  in  obtaining  the 
velocity  components  is  less  than  two  percent.  For  the  five-hole  yaw-probe, 
measurements  of  pressure  coefficients,  however,  appreciable  errors  were 
found.  The  test  result  showed  that,  for  the  pitch  and  yaw  angles  between  zero 
and  15  degrees,  the  pressure  coefficients  have  values  ranging  from  0.006  up  to 
0.02,  in  contrast  to  the  actual  value  of  zero.  These  variations,  at  the  first 
glance,  may  seem  small  enough,  but  they  can  have  a  significant  effect,  especi¬ 
ally  near  the  midship,  where  the  actual  pressure  coefficients  are  small  and  of 
the  same  order  as  the  numerical  errors.  For  this  reason  it  is  recommended 
that  the  5-hole  probe  measurements  of  pressure  coefficient,  especially  around 
the  midship,  be  disregarded. 

From  the  discussion  in  this  section,  it  may  be  concluded  that  the  overall 
physical  and  numerical  errors,  except  for  five-hole  probe  measurements  of 
pressure  coefficients,  are  small,  and  the  present  measurements  are  of  accept¬ 
able  accuracy  within  an  uncertainty  of  a  few  percent.  It  should  be  noted  that 
none  of  the  data  was  corrected  for  the  aforementioned  physical  or  numerical 
errors. 
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Figure  10.  Comparison  of  Ct's  vs.  Froude  Nuqiber.  (Free  Condition) 
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2.  Total-Resistance  Measurements 


Total-resistance  coefficients  C^,  obtained  from  measurements,  are  given 
in  Tables  1  and  2  from  tests  at  fixed  and  f ree-conditions,  respectively. 
These  results  are  compared  in  Figures  10  and  11  with  the  other  experimental 
results;  those  for  a  2.5-meter  model  tested  at  the  University  of  Tokyo  (UT),  a 
4-m  model  tested  at  the  Ship  Research  Intitute  of  Japan  (SRI),  and  a  6-m  model 
tested  at  the  Ishikawaj ima-Harima  Heavy  Industries  Co.  Ltd.  (IHI).  The  latter 
model  has  been  tested  only  in  the  free  condition.  To  have  a  meaningful  com¬ 
parison,  all  the  results  were  corrected  to  a  standard  temperature  of  18.3°C 
for  the  fixed-condition  tests  and  13. 3° C  for  the  f ree-condition  tests.  These 
results  show  a  similar  oscillatory  variation  with  Froude  number,  i.e.  a  signi¬ 
ficant  hump  and  hollow  in  the  range  0.31  <  Fr  <  0.32  and  0.34  <  Fr  <  0.35, 
respectively.  There  are  some  additional  humps  at  Froude  numbers  of  0.21  and 
0.25,  and  hollows  at  Froude  numbers  of  0.22  and  0.27. 

The  comparison  of  results  between  free  and  fixed  conditions  shows  that 
the  humps  and  hollows  for  both  cases  occur  at  the  same  Froude  numbers,  but  the 
values  of  Ct  with  the  free  condition  are  greater  than  those  with  the  fixed 
condition,  by  one  to  four  percent  for  Froude  numbers  less  than  0.28,  and  4  to 
10  percent  for  Froude  numbers  greater  than  0.28. 

The  residuary  resistance  Cw  for  the  fixed  condition,  from  the  three 
towing  tanks,  also  have  been  shown  in  Figure  10.  These  values  have  been 
calculated  using  the  form-factor  formula 

Ct  =  (1+k)  Cf  +  Cw  (8) 

where  is  the  frictional-resistance  coefficient  for  a  flat  plate,  derived 
from  the  Schoenherr  (1932)  formula,  and  k  is  a  form  factor,  which  has  the 
values  of  k  =  0.10  for  the  Iowa,  k  =  0.065  for  SRI,  and  k  =  0.05  for  'JT. 
These  values  were  selected  so  as  to  obtain  the  best  agreement  between  the  Cw's 
of  the  three  tanks.  The  higher  value  for  the  Iowa  towing  tank  is  due  to  the 
fact  that  this  tank  has  the  smallest  width  among  the  tanks,  and  as  is  shown  by 
Landweber  and  Nakayama  (1975),  the  measured  values  of  total-resistance  coef¬ 
ficients  are  higher  than  the  measured  values  in  a  wider  tank  for  Fr  >  0.30  due 
to  blockage  effect  of  the  towing-tank  walls.  The  form-factor  procedure  gives 
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a  good  estimate  of  the  viscous  resistance  for  the  Wigley  hull  due  to  its  low 
block  coefficient.  This  has  been  shown  by  Ju  (1983),  who  showed  that  the 
viscous-resistance  coefficient  obtained  from  a  wake  survey  is  weakly  dependent 
upon  the  Froude  number  for  the  Wigley  hull. 

3.  Pressure-Distribution  Mesurements 

The  results  of  the  pressure-distribution  measurements  on  the  Wigley  hull 

along  waterlines  j Z/H |  =  0.2  and  0.6  are  shown  in  Table  3*  The  pressure 

coefficient  Cp  is  defined  as  Cp  =  (Pg  -  PJ/  ^  puf ,  where  Pg  is  the  pressure 

on  the  hull  and  is  the  speed  of  the  model.  The  experimental  procedure  is 

carried  out  in  such  a  way  that  P  =  0  in  the  above  relation.  These  results 

00 

are  compared  with  the  potential-flow  calculations  by  Guilloton's  (1951)  method 
and  are  shown  in  Figures  12  through  15.  These  comparison  show  fair  agreement 
between  calculations  and  measured  values  over  most  of  the  hull  except  near  the 
bow  and  stern.  This  may  be  partly  due  to  a  logarithmic  infinity  in  the 

Guilloton  calculations  at  the  bow  and  stern  for  negative  values  of  Z/H,  i.e. 
beneath  the  free  surface.  At  the  free  surface,  the  logarithmic  part  of  the 

calculated  functions  vanishes,  so  that  the  measured  wave  profile  along  the 

hull  is  in  very  good  agreement  with  the  result  of  calculations  there.  An 

important  feature  of  these  figures  is  the  fact  that  the  pressure  distribution 
on  the  hull  closely  follows  the  trend  of  the  wave  profile  at  the  free  surface, 
i.e.,  shows  the  same  hollows  and  humps.  This  shows  that  the  Froude  number  has 
an  important  effect  on  the  boundary  layer  and  its  development  on  ship  hulls. 

4.  Velocity  Distribution 

4.1.  General 

A  sample  of  tne  output-data  file  for  five -hole -probe  measurements  is 
shown  in  Table  4.  There  are  a  few  points  about  this  output  file  which  should 
be  mentioned  in  order  to  have  a  better  grasp  of  these  measurements.  The  first 
column  shows  the  distance  of  measured  points  from  the  hull  along  the  y-axis, 
i.e.,  at  constant  waterline  planes.  The  first  point  of  these  measurements, 
where  the  5-hole  probe  is  in  contact  with  the  hull,  is  denoted  as  the  zero- 
distance  point,  and  the  rest  of  the  points  are  referenced  relative  to  this 
point.  The  actual  distance  of  the  probe  axis  from  the  hull  at  this  point  is 
2«15  mm.  Columns  5,  6,  and  7  show  the  nondimensional  velocity  components  u. 


33 


II 

I! 


il 


_ i  1 1  -^r*  in  OJ  *r  CO  'CO  C*~»  OJ  O'  CM  CD 

•=»  il  — ■*  -*r  C\j  -z=.  ru  ?J  — «  <r  r-o— « cn  =a  -o 

'i~  » i  ->o  -o  -r  co  co  m  ?•'.  co  co  r-".  a  < — i  •— *  c=j  — • 

H  rza  *=•  ro  co  '=S**^»  *ro  «•-:=»  cs*  c=.cor=o  co  co 

co  U 

iio^-zz.cur^w?.  —  “ 


I 


t  i 


a  ""*"  I  -,r*  -T-  ,TD 
•«»•  J'OszaN  M 
*•  *0  ’•  j<=>'ca'-0  -r 


* I  zzj  o  co  o  «ro  ^  >— • 


Qu 

O 


CsS  II  -TZ3  It 


*  -C=>  uCi  .  rj  C  ~M 

*  *D  ‘U.-I  c  -»  oj 

"-Ti^OOo 


Cid 

i — 

CO 


H  «j*  <>“•  LO  «"S.  Cl'-  <T 

1 1  *7-  Lr  *  75  lo  -  o  O'-  r  >v 
ii  r--  co  .  o  •— *  ->-•  '•.'■•  cj 


•>  Ufl  O'-  rzz i  O'*-  CO  U"1  O*. 
Ci"-  r-  C^-  c-'i  -rr  •— 1  c*’i 
T*  0>  -O  -i*  r-~ 


1 1  co  *co  co  co  *T«  co  co  o  era  co  •■ 


erv  co  -r-a  co  —  •--  co  co  era  — •  rro 
^  lioooooooc-i  <zz-  u*i 


03  £  o 
C*-  OO  >-  < 


sCJnCJ  vfO>siOri  - 

a  — «  3  0[v3<r>-H; 

a  .iO  cj  cj  o:0'0ro  oj  ■« 


i  co  co  co  co 


<C 

CO 


co  il  -q-rv<]C3ij-ii>.,J-ii>.^c3un-rv0^r's. 

i_n  1 1  x  o-  co  ro  oj  oj  o-  oj  o  oj  ~-0  — ■— *  m 
on  u  o j  =3  <=>  o  un  ~0  un  rO  o  m  oj  '==*  co  <^j  r *d 

.  II  ^,'M^hOO'^^'^OCoO’OOO 

co  II  .  •  •  • 

II  oO  ra  co  co  co  co  co  <o<o  co  co  era  co  c=> 


CD 


m  H  0'0c^-ocs.o0'0cooo'«oococo 

— i  1 1  -r— <  — <  lh  O-  "-Q  — *  "O  -o  O  **■»♦  O  -o  ■<*"  -r- « 

,-D  1 1  ^rOrrao2<rO'M^r<i'„-i>.un'r 


-/3  II  03  CO  CD  O  O  o  o  O  co  O  O*.  Oca  O'-  C>- 

-JO  II  O^OC^OJcoOO'~O^T*'Or-DO'JD'-OC>- 

oj  m  un-^r  *r>a-jD‘oo'»^oco<T-^o^r<i 


oc 

DD 

CO 


_ I  II  LnoLnooooi 

\  ii  o-'OCQ'S  — o -<r o j < 
X  II  ••.... 
OJ  1 1  co  CO  CO  co  — o  CO  CO  ^ 


34 


Table  3.  Measured  pressure-distribution  coefficient 


Cm.  OF  WATER 


2X/L 

Figure  12.  Comparison  of  measured  and  calculated  pressure 
distribution  (Fr  =  0.266). 
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Cm.  OF  WATER 


2X/L 

Figure  13.  Comparison  of  measured  and  calculated  pressure 
distribution  (Fr  =  0.313). 
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OF  WATER 


2X/L 

figure  14.  Comparison  of  measured  and  calculated  pressure 
distribution  (Fr  =  0.350). 
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Cm.  OF  WATER 


Figure  15.  Comparison  of  measured  and  calculated  pressure 
distribution  (Fr  =  0.400). 
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BOUNDARY  LAYER  MEASUREMENTS  FOR  UIGLEY  HULL 


LENG7H=304.3  CM  BEAH=30 . 48  CM  DRAF7=19.05  CM 
FROUDE  NO. =0.400  2X/L=0.866  Z7H=0.4 


DIST  (nn) 

H/H  ret 

Cp  Q/Q  ref 

0.000 

0.269 

0.040 

0.229 

2.005 

0.325 

0.020 

0.305 

4.010 

0.382  • 

-0.001 

0.332 

6.015 

0.421  ■ 

-0.005 

0.427 

8.020 

0.463 

-0.014 

0.477 

10.025 

0.497 

-0.022 

0.519 

12.030 

0.529 

-0.023 

0.552 

14.026 

0.557 

-0.026 

0.583 

16.041 

0.577 

-0.019 

0.596 

13.046 

0.655 

-0.020 

0.625 

■20.051 

0.629 

-0.015 

0.644 

•22.056 

0.652 

-0.017 

0.669 

26.066 

0.682 

-0.015 

0.697 

30.076 

0.720 

-0.021 

0.741- 

34.096 

0.772 

-0,025 

0.797 

38.096 

0.827 

-0.024 

0.352 

42.107 

0.375 

-0.Q22 

0.897 

46.117 

0.919 

-0.018 

0.937 

50.127 

0.944 

-0.016 

0.960 

60.152 

0.986 

-0.010 

0.996 

70.178 

1.005 

-0.003 

1.014 

80.203 

1.012 

-0.007 

1.019 

90.223 

1.311 

-0.008 

1.019 

110.279 

1.003 

-0.011 

1.014 

130.330 

0.999 

-0.015 

1.014 

y/Q  ref 

v/Q  ref 

u/Q  ref  TKETA 

PSI 

0.222 

0.026 

0.043 

6.3  ■ 

-12:2 

0.295 

0.045 

0.062 

8.4 

-11.9 

0.369 

0.060 

0.077 

9.3 

-11.3 

0.413 

0.065 

0.085 

9.1 

-11.7 

0.463 

0.069 

0.095 

8.4 

-11.6 

0.505 

0.072 

O'.  101 

8.1 

-11.3 

0.537 

0.072 

0.107 

'7.7 

-11.2 

0.568 

0.074 

0.111 

7.5 

-11.0 

0.580 

0.074 

0.114 

7.3 

-11. 1 

0.609 

0.078 

0.117 

7.3 

-10.9 

0.623 

0.079 

0.119 

7.2 

-10.8 

0.653 

0.031 

0.121 

7.1 

-10.5 

0.631 

0.081 

0.122 

6.3 

-10.2 

0.726 

0.031 

0.125 

6.4 

-9.3 

0.783 

0.082 

0.126 

6.0 

-9,2 

0.938 

0.085 

0.126 

5.3 

-9.6 

0.884 

O.0S8 

0.125 

5.7 

-8.1 

0.925 

0.089 

0.123 

5.5 

-7.6 

0.950 

0.090 

0.121 

5.4 

-7.3 

0.989 

0.086 

0.114 

5.0 

-6.6 

1.008 

'0.082 

0.109 

4.7 

-6.2 

1.011 

0.077 

0.104 

4.4 

-5.8 

1.011 

0.074 

0.099 

4.1 

-5.5 

1.008 

0.071 

0.091 

3.9 

-5.0 

1.009 

0.067 

0.034 

3.5 

-4.5 

Table  4.  Sample  of  output  file  on  tape. 
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V,  and  w  along  the  x,  y,  and  z  axes,  respectively,  u  >  o  denotes  the  velocity 

direction  towards  the  stem;  w  >  o,  represents  the  velocity  direction  from 
keel  toward  free  surface;  and  v  >  o  denotes  the  flow  direction  toward  tne 

hull.  Column  4  shows  the  nondimens ional  magnitude  of  the  velocity  vector,  and 
column  3  gives  the  pressure  coefficients  measured  by  the  five-hole  probe. 
Column  2  represents  the  total  head  which  is  obtained  from  columns  3  and  4. 

Columns  8  and  9  show  the  pitch  and  yaw  angles,  respectively,  which  are  defined 

in  Figure  2.  The  complete  set  of  data  for  all  four  Froude  numbers,  is  ar¬ 
chived  on  magnetic  tape  and  is  available  from  the  IIHR  tape  library. 

4.2.  Velocity  Distribution  on  the  Hull 

Figures  16  through  20  show  typical  distributions  of  the  u-component  of 
the  mean  velocity  vector.  It  is  seen  that  the  flow  is  of  a  mild  adverse- 
pressure-gradient  type.  The  profiles  also  seem  to  have  a  characteristic 

concavity  near  the  wall.  This  is  due  to  the  beginning  of  the  wake-like  behav¬ 
ior  of  the  outer  (fully  turbulent)  layer.  These  figures  show  that  the 

boundary  layer  is  rather  thin  over  most  of  the  hull,  except  near  the  stem 
where  the  thickness  becomes  of  the  order  of  four  percent  of  the  model 

length  (f  =  0.04);  but  still  it  may  be  regarded  as  a  ’thin'  boundary  layer. 
There  are  no  drastic  changes  in  the  profile  in  the  x-direction.  Also,  there 
is  no  reversal  of  the  u-component  velocity  even  at  the  stem,  which  indicates 
that  the  flow  does  not  separate  over  the  Wigley  hull.  An  interesting  feature 
of  the  boundary  layer  which  can  be  seen  from  these  figures  is  its  rate  of 
growth.  The  boundary-layer  thickness  at  the  waterline  j Z/H |  -  0.2  grows  along 
the  hull;  see  Figure  16;  however,  the  rate  of  growth  decreases  at  the  deeper 
waterlines,  as  is  seen  at  |z/h|  =  0.4  (Figure  17)  and  becomes  almost  constant 
at  | Z/H |  =  0.6;  see  Figure  18.  Near  the  keel,  the  rate  of  growth  reverses,  as 
is  shown  in  Figures  19  and  20  for  | Z/H |  =0.8  and  | Z/H |  =  1.0,  respectively. 
There  the  boundary-layer  thickness  is  largest  amidships  and  diminishes  along 
the  hull  until  the  stem  where  it  becomes  almost  zero.  This  phenomenon  can  be 
attributed  mainly  to  the  geometry  of  the  Wigley  hull.  The  streamlines  from 
bow  to  midship  converge  toward  the  keel,  causing  the  boundary  layer  to  become 
thick  near  the  keel  at  midship.  From  there  to  the  stem,  the  reverse  of  the 
aforementioned  phenomenon  occurs;  i.e.  the  streamlines  diverge  from  the  keel, 
causing  the  boundary  layer  to  become  thinner  as  the  flow  approaches  the  stem. 
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Consider  now  the  distribution  of  the  v-  and  w-eomponents  of  the  mean 
velocity.  The  w-component  at  all  the  sections  is  positive,  i.e.,  its  direc¬ 
tion  is  from  the  keel  towards  the  free  surface.  This  can  be  explained  by 
noting  that  the  wave  profiles  on  the  hull  have  a  trough  at  about  the  raid- 
ship.  As  will  be  shown  in  Chapter  IV,  the  streamlines  on  the  hull  closely 
follow  the  trend  of  the  surface-wave  profiles,  so  that  the  flow  from  midship 
to  the  stern  is  also  upward;  i.e.,  the  w  component  will  be  positive  in  that 
region. 

The  direction  of  the  v-component  of  the  velocity,  for  the  most  part,  is 
toward  the  hull.  This  should  be  expected  since  the  flow  from  midship,  where 
the  hull  is  thickest,  converges  on  the  both  sides  of  the  model  toward  the 
stern.  This  implies  that  the  velocity  vector  has  a  component  in  the  negative 
direction  of  the  y-axis  on  the  starboard  side  where  the  measurements  were 
taken.  Figures  21  through  28  show  some  typical  vector  plots  of  the  v-w  com¬ 
ponents.  An  interesting  feature  of  these  plots  is  shown  at  the  |  Z/H |  =  0.2 
waterlines  which  are  very  close  to  the  free  surface.  Since  —  =  at  the  free 
surface,  the  slope  of  the  velocity  vector  at  this  waterline  indicates  the 
effect  of  the  boundary  layer  on  the  lateral  slope  of  the  free  surface. 

5.  Wake  Measurements 

The  mean-velocity  components  were  measured  at  four  stations  from  the 
stern  to  30.48  cm.  Figures  29  through  32  show  the  typical  distribution  of  the 
u  components  of  the  velocity  vector.  The  profiles  are  seen  to  be  symmetric. 
The  point  of  the  minimum  velocity  occurs  at  the  wake  centerline.  Since  the 
measurements  were  not  carried  out  far  enough,  it  was  not  possible  to  investi¬ 
gate  the  decay  of  three  dimensionality,  or  the  asymptotic  behavior  of  the 

wake.  Figures  33  through  40  show  the  distribution  of  the  v-w  components 

2x  2x  1 
at  —  =  0  and  —  =  1.2,  for  various  Froude  numbers. 

6.  Effects  of  Froude  Number  on  Boundary  Layer 

The  practical  importance  of  free-surface  phenomena  in  viscous  flow  around 
ships  has  been  reviewed  by  Baba  (1981).  Also  Landweber  (1971),  and  Landweber 
and  Patel  (1979)  have  discussed  in  detail  the  various  effects  of  the  free 
surface,  such  as  the  free-surface  boundary  layer,  bilge  vortices,  breaking  of 
the  bow  wave  and  separation  of  free  surface.  In  this  section,  we  shall  point 
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out  the  variation  of  the  measurement  data  caused  by  changes  of  the  Froude 
number,  and  also  compare  the  present  results  with  the  measurements  obtained  by 
Sarda  (1985)  on  a  double  model  of  the  Wigley  hull  of  the  same  dimensions,  in 
the  IIHR  wind  tunnel. 


The  first  characteristic  showing  the  influence  of  the  Froude  number  is 
the  pressure  distribution  over  the  hull.  In  the  case  of  zero  Froude  number, 
the  measured  pressure  coefficient  is  a  nearly-symmetric  curve  about  the  mid¬ 
ship  at  any  waterline  plane.  In  the  case  of  nonzero  Froude  numbers,  however, 
the  measured  pressure  coefficients  are  strongly  affected  by  the  wave  profile 
along  the  hull,  and  closely  follow  the  trend  of  the  free  surface  as  is  shown 
in  Figures  12  through  15.  This  fact  is  re-emphasized  in  Figures  41  through 
44.  These  figures  show  th^  comparison  of  the  measured  pressure  coefficient 

and  wave  profile  with  1  -  (jr^ — )  ,  where  Qg  is  the  magnitude  of  the  velocity 

ref 

vector  at  the  edge  of  the  boundary  layer.  This  consequence  of  the  constant 
free-surface  pressure  at  nonzero  Froude  numbers  significantly  affects  many 
characteristics  of  the  boundary  layer,  such  as  the  rate  of  growth,  the  magni¬ 
tude  of  the  cross  flow,  the  piezometric  pressure  gradients,  and  the  velocity 
at  the  edge  of  the  boundary  layer.  Figure  45  shows  the  comparison  of  the 
displacement  thickness  along  the  hull  at  | z/H |  =  0.2  for  different  Froude 
numbers.  There  is  a  combination  of  Reynolds-and-Froude-number  effects  on  the 
rate  of  growth  of  the  boundary  layer  along  the  hull.  The  Reynolds  number 
influences  the  boundary  layer  through  viscosity  whereas  the  Froude  number 
imposes  itself  mainly  through  pressure  distribution.  At  the  first  station, 
2x/L  =  0.0,  the  ratio  of  displacement  thickness  at  Fr  =  0.266  to  the  one  at  Fr 
=  0.400  is  about  2.  This  ratio  for  comparable  Reynolds  numbers  of  a  double 
model  is  about  1.2.  This  is  the  first  indication  of  Froude-number  effect  on 
boundary -layer  development.  The  rate  of  growth  of  boundary  layer  for 
different  Froude  numbers  then  exactly  follows  the  trend  of  the  pressure  dis¬ 
tribution  at  z/H  =  0.2  (Figures  12  through  15,  or  Table  3).  From  midship  to 
2x/L  =  0.4,  the  flow  is  decelerating  for  all  the  Froude  numbers  except  Fr  = 
0.266.  Therefore  the  displacement  thickness  has  the  smallest  rate  of  growth 
for  Fr  =  0.266  in  that  region.  From  2x/L  =  0.4  to  2x/L  -  0.6,  displacement 
thickness  has  its  largest  value  for  Fr  =  0.350.  This  is  because  flow  over  the 
Wigley  hull  at  this  station  has  the  greater  rate  of  deceleration  in  comparison 
with  the  other  Froude  numbers,  and  therefore  the  boundary  layer  becomes  thick 
faster.  Between  2x/L  =  0.6  and  2x/L  =  0.8  stations,  only  flow  at  Fr  =  0.313 
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is  accelerating  and  the  displacement  thickness  at  this  Froude  number  has  a 
smaller  value  -in  comparsion  with  the  other  Froude  numbers.  Finally,  near  the 
stern,  the  displacement  thickness  at  Fr  =  0.266  has  the  largest  rate  of  growth 
because  it  has  the  largest  unfavorable  pressure  gradient  and  flow  decelerates 
at  a  faster  rate  than  at  the  other  Froude  numbers.  Figures  25  through  28  show 

the  vector  plot  of  crossflow  at  2x/L  =  0.866.  The  wave  profile  on  the  hull  at 

this  station  has  slopes  of  0.051,  0.034,  0.032,  and  0.046  for  Froude  numbers 

of  0.266,  O.313,  n.350,  and  0.400,  respectively.  Comparison  of  these  figures 

shows  that  the  magnitude  of  the  crossflow  is  directly  related  to  the  slope  of 
the  free  surface.  For  Froude  numbers  of  0.313  and  0.350,  where  the  slope  is 
almost  the  same,  the  crossflow  vectors  are  very  similar  and  their  magnitudes 
are  smaller  than  those  at  Froude  numbers  of  0.266  and  0.400.  As  a  matter  of 
fact,  at  Fr  =  0.266,  the  crossflow  is  larger  than  at  the  other  Froude 
numbers.  Figures  46  through  49  show  the  comparison  of  the  yaw  angle  V ,  which 
is  directly  proportional  to  the  magnitude  of  the  crossflow  component  'of  the 
velocity  vector,  at  two  stations  2x/L  =  0.8  and  0.933,  for  zero  and  nonzero 
Froude  numbers.  The  comparison  shows  that  the  yaw  angles  for  nonzero  Froude 
numbers  are  at  least  200  percent  of  the  comparable  angles  for  the  double 
model.  This  comparison  also  show  that  the  effect  of  Froude  number  on  the 
boundary  layer  is  felt  primarily  near  the  free  surface,  as  is  shown  in  Figures 
46  and  48  for  | z/H |  =  0.2.  Figures  47  and  49  show  the  yaw-angle  distribution 
across  the  boundary  layer  for  |z/H|  =  0.8.  Note  that  the  yaw  angles  for 
different  Froude  numbers  almost  collapse  into  a  single  curve  near  the  keel, 
but  still  they  are  almost  twice  the  comparable  angles  for  the  double  model. 
The  effect  of  the  free  surface  on  cross-flow  extends  into  the  wake  region  of 
the  model.  Figures  37  through  40  show  the  cross-flow  plots  at  2x/L  =  1.2  for 
different  Froude  numbers.  Note  the  distinct  S-shaped  profile  which  occurs  at 

Froude  numbers  0.266,  0.313,  and  0.350,  while,  at  Fr  =  0.400  the  cross-flow 
does  not  show  such  a  behavior* 


43 


4^ 


140. 


FR-0.  266 
Z/H-0.  2 


E 

E 

D 

X 

W 

X 

h* 


120. 


100.  - 


80.  - 


o 

o 

II 

\ 

X 

04 


04 

O 

II 

A 

\ 

X 

04 


O 

II 

\ 

X  • 

04 


n 

o 

II 

\ 

X 

04 


« 

00 

866 

* 

• 

• 

• 

• 

o 

o 

o  • 

O 

II 

II 

II 

11 

1 

\ 

\ 

# 

\  • 

X  * 

X 

X 

X 

04 

04 

04 

04 

• 

« 

« 

o 

q: 

U- 

LlI 

u 

< 

H 

0) 

1-4 

Q 


40. 


Figure  16.  Measured  streamwise  velocity  component  (|Z/H|  -  0.2). 
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Figure  17.  Measured  streamwise  velocity  component  ( J Z/H |  =  0.4). 
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Figure  18.  Measured  streamwise  velocity  component  ( | Z/H |  -  0.6). 
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Figure  19.  Measured  streamwise  velocity  component  ( | Z/H |  -  0.8). 
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Figure  20.  Measured  streamwise  velocity  component  (]Z/H|  =  1.0). 
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Figure  21.  Vector  plot  of  w-v  components  ( 2 x/L  -  0.5,  Fr  -  0.266). 


Figure  22.  Vector  plot  of  w-v  components  (2x/L  =  0.5,  Fr  _  0.313) 
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Figure  23.  Vector  plot  of  w-v  components  (2x/L  =  0.5,  Fr  =  0.350). 
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Figure  24.  Vector  plot  of  w-v  components  (2x/L  =  0.5,  Fr  =  0.400) 
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Figure  25.  Vector  plot  of  w-v  components  (2x/L  =  0.86,  Fr  =  0.266). 
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Figure  26.  Vector  plot  of  w-v  components  (2x/L  =  0.86,  Fr  =  0.313). 
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Figure  27.  Vector  plot  of  w-v  components  (2x/L  =  0.86,  Fr  =  0.350). 
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Figure  28.  Vector  plot  of  w-v  components  (2x/L  =  0.86,  Fr  =  0.400). 
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Figure  29.  Measured  streamwise  velocity  component  in  the  wake 
(  z/H  =  0.2,  Fr  =  0.350. 
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Figure  30.  Measured  streamwise  velocity  component  in  the  wake 
4  ( | z/H  =  0.4,  Fr  =  0.400). 
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Figure  31.  Measured  streamwise  velocity  component  in  a  cross 
section  (2x/L  =  1.1,  Fr  =  0.266). 
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Figure  32.  Measured  streamwise  velocity  component  in  a  cross 
section  (2x/L  =  1.1,  Fr  =  0.313). 
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Figure  33.  Vector  plot  of  w-v  component  (2x/l  =  1.0,  Fr  =  0.266). 
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Figure  34.  Vector  plot  of  w-v  component  (2x/L  =  1.0,  Fr  =  0.313). 
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Figure  35.  Vector  plot  of  w-v  component  (2x/L  =-1.0,  Fr  =  0.350) 
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Figure  36.  Vector  plot  of  w-v  component  (2x/L  *  1.0,  Fr  =  0.400). 
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Figure  37.  Vector  plot  of  w-v  component  (2x/L  =  1.2,  Fr  =  0.266) 
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Figure  33.  Vector  plot  of  w-v  component  (2x/L  =  1.2,  Fr  =  0.313) 
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Figure  39.  Vector  plot  of  w-v  component  (2x/L  =  1.2,  Fr  =  0.350). 


67 


Z/H 

-0.8  -0.6  -0.4  "0* 2 


o 


DISTANCE  FROM  THE  WAKE  CENTERLINE  (mm) 

Figure  40.  Vector  plot  of  w-v  component  (2x/L  =  1.2,  Fr  =  0.400). 
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Figure  44.  Comparison  of  measured  Cn  and  1  -  (CL/Qnaf)^ 
(Fr  =  0.400).  P  e  rer 
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Figure  45.  Comparison  of  displacement  thickness  along  the  hull 
( | z/H |  =  0.2). 
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Figure  47.  Comparison  of  the  yaw  angles  across  the  boundary  layer 
(2x/L  =  0.8,  I z/H I  =  0.8). 
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Figure  49.  Comparison  of  the  yaw  angles  across  the  boundary  layer 
(2x/L  =  0.933,  | z/H |  =  0.8). 
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CHAPTER  IV 

BOUNDARY-LAYER  CALCULATIONS  FOR  THE  WIGLEY  HULL 


1.  General 

The  problem  of  predicting  ship  viscous  resistance  has  intrigued  physi¬ 
cists,  mathematicians,  and  engineers  for  many  years.  The  inadequacy  of  the 
so-called  Froude  method,  in  which  the  viscous  resistance  is  taken  as  the  fric¬ 
tional  resistance  of  an  equivalent  flat  plate  with  an  empirical-shape  factor 
correction,  has  been  demonstrated  experimentally  (Landweber  and  Wu,  1963). 
Reliable  and  accurate  methods  for  determining  the  viscous  resistance  should 
necessarily  involve  the  solution  of  the  three-dimensional  boundary-layer 
equations.  In  view  of  the  difficulties  of  calculating  the  pressure  distribu¬ 
tion  on  a  ship  form  at  a  nonzero  Froude  number,  however,  the  computation  of 
the  boundary-layer  development  on  ship  hulls  seems  a  formidable  task.  Due  to 
this  complexity,  three  dimensional  boundary  layers  on  ship  hulls  are  often 
calculated  by  considering  the  double  model,  i.e.,  the  zero-Froude-number 
case.  The  influence  of  the  free  surface  on  the  boundary  layer  is  therefore 
neglected.  A  review  of  these  methods  is  given  by  Landweber  (1971),  and 
Landweber  and  Patel  (1979).  The  aforementioned  simplification  suffers  from 
several  drawbacks.  First,  the  pressure  distribution  on  the  double  body  is 
independent  of  Froude  number.  In  the  nonzero  case,  however,  the  pressure 
distribution  is  a  function  of  the  Froude  number.  Consequently,  the  boundary 
layer  over  the  hull  is  not  modeled  correctly  by  double-body  calculations. 
Secondly,  the  changes  in  the  wetted  area  of  the  hull,  associated  with  differ¬ 
ent  wave  patterns,  are  ignored  and  therefore  the  influence  of  Froude  number  on 
viscous  resistance  cannot  be  calculated.  Finally,  any  local  separation  due  to 
wave  steepness  along  the  hull  (Chow,  1967),  cannot  be  predicted. 

As  a  first  step  toward  taking  a  wavy  free  surface  into  account, 
Shahshahan  and  Patel  (1983)  considered  the  boundary  layer  along  the  intersec¬ 
tion  between  the  hull  and  the  free  surface  and  demonstrated  the  phenomenon  of 
local  wave-induced  separation  at  the  free  surface.  In  this  chapter,  their 
method  has  been  employed  to  calculate  the  development  of  the  boundary  layer  on 
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the  Wigley  hull.  The  numerical  results  have  then  been  compared  with  the 
experimental  results  obtained  in  the  towing  tank. 

2.  Coordinate  System  and  Governing  Equations 
2.1.  Coordinate  System 

Consider  the  ship  hull  with  its  surface  is  given  by  y  =  ±  f(x,z)  in 
rectangular  Cartesian  coordinates  (x,y,z),  with  increasing  values  of  x  in  the 
direction  opposite  to  the  ship's  forward  motion,  z  vertically  upward,  and  the 
origin  at  midship  at  the  undisturbed  free  surface.  Denote  any  streamline  on 
the  hull  by  z  =  g(x).  This  streamline  is  assumed  to  be  known  either  from 
experiment  or  inviscid-flow  theory.  For  the  hull  boundary  layer  along  this 
streamline  an  orthogonal  coordinate  system  (5  ,n  ,£  )  has  been  chosen,  in  which 
£  is  measured  along  the  streamline,  n  is  normal  to  the  hull  (n  =  o  at  the 
hull  surface),  and  £  is  the  third  member  of  this  right-handed,  orthogonal 
coordinate  system.  See  Figure  50.  It  should  be  noted  that  this  coordinate 
system  is  triply  orthogonal  only  on  the  hull  surface,  and  it  is  not  exactly 
orthogonal  off  the  hull  since  the  £  and  £  lines  do  not  coincide  with  the  lines 
of  principal  curvature  of  the  surface.  The  departures  from  orthogonality  of 

the  coordinate  system,  however,  will  be  small  since  the  boundary  layer  is 
assumed  to  be  thin. 

In  the  present  calculation,  the  measured  wave  profile  at  the  hull  was 
used  as  the  first  streamline.  Other  streamlines  along  the  hull  below  the  free 
were  calculated  by  using  Guilloton's  (1948,  1951)  method.  This  method  has 
been  used  extensively  by  many  researchers,  e.g.,  Emerson  (1967),  Webster  and 
Huang  (1970),  Gadd  (1973),  and  Adachi  and  Hinatsu  (1981),  and  indications  are 
that  the  resulting  streamlines  are  in  good  agreement  with  measurements.  The 
single  drawback  is  that  only  a  few  streamlines  can  be  obtained  by  this  proced¬ 
ure.  Figures  51  through  54  show  the  measured  wave  profile  and  four  stream¬ 
lines  which  were  obtained  for  the  Wigley  hull  at  Froude  numbers  0.266,  0.313, 
0.350,  and  0.400,  respectively. 
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Figure  50.  ship  (x.y.z)  and  bounda 


ry  layer  (f,n,c)  coordinates. 
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Figure  53  Measured  wave  profile  and  Gullloton's  streaml Ine  (Fr  ®  0.350). 
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Figure  54  Measured  wave  profile  and  Guilloton’s  streaml Ine  (Fr  B  0.400) 


2.2.  Governing  Equations 

One  of  the  most  common  assumptions,  which  has  been  used  extensively  in 
boundary-layer  calculations  applicable  to  ship  forms,  is  the  assumption  of 
small  crossflow,  i.e.,  the  direction  of  flow  within  the  boundary  layer  devi¬ 
ates  by  only  a  small  angle  from  the  direction  of  the  streamline  at  the  outer 
edge  of  the  boundary  layer.  The  validity  of  this  assumption,  however,  has 
been  questioned;  Landweber  (1971).  The  combination  of  the  hull  geometry  and 
the  secondary  flow  generated  by  the  wave  profile  along  the  hull,  particularly 
around  the  bow  and  stem  where  the  wave  profile  is  the  steepest,  results  in  a 
large  crossflow  velocity  component  for  most  hulls.  Therefore  the  small- 
crossflow  assumption  seems  unsuitable.  For  the  Wigley  hull,  however,  this  is 
not  the  case.  The  boundary-layer  measurements  obtained  for  this  model  show 
that  the  crossflow  angle  (6  =  tan  reaches  a  maximum  value  of  about  eight 
degrees  for  the  steepest  wave  profile.  Thus  one  may  assume  that  the  boundary 
layer  over  the  Wigley  hull  can  be  calculated  by  using  the  small-crossflow 
equations  (see  Nash  and  Patel  (1972)) 
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puv  and  pvw  are  Reynolds  stresses,  v  is  the  Kinematic  viscosity,  and  U0  is 
the  velocity  in  the  5  direction  outside  the  boundary  layer. 

The  metric  coefficients  h^ ,  h  and  the  coordinate  curvatures 
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can  be  related  to  the  hull  geometry  f(x,z),  and  the  streamline  equation  g(x) 
using  "ths  moHiod  of*  Miloh.  nnd  Pulo!  (1973)  ns  follows j 
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Equation  (17)  expresses  the  orthogonality  of  the  ?  and  K  directions. 


The  assumptions  of  small  crossflow  enable  the  boundary  layer  along  each 
streamline  to  be  calculated  independently.  This  argument  is  also  valid  at  the 
free  surface.  The  assumption  that  the  deviations  of  the  t,  coordinate  line 
from  the  normal  to  the  free  surface  across  the  boundary  layer  are  small,  i.e., 
the  deflections  of  the  free  surface  across  the  boundary  layer  are  small,  then 
leads  to  the  usual  small-crossflow  approximations  and  enables  the  boundary 
layer  along  the  wave  profile  to  be  calculated  independently  of  the  flow  below 
the  free  surface. 


3.  Solution  Procedure  and  the  Results 
3.1.  Solution  Procedure 

The  small-crossflow  assumption  uncouples  the  continuity  equation  (10)  and 
the  £ -momentum  equation  (9)  from  the  ? -momentum  equation  (11).  Hence  equa¬ 
tions  (9)  and  (10)  can  be  solved  for  U  and  V,  and  then  the  crossflow  W  can  be 
determined  from  equation  (11).  These  equations  have  been  solved  by  the  Crank— 
Nieolson  method  of  Chang  and  Patel  (1975).  The  Reynolds  stresses  are  obtained 
from  the  two-layer  eddy -viscosity  model,  and  transition,  wherever  applicable, 
is  simulated  by  an  intermittency  function.  See  Cebeci  and  Smith  (1976). 

Boundary-layer  calculations  along  the  wave  profiles  on  the  hull  were 
started  with  the  Blasius  solution  a  short  distance  downstream  of  the  bow  and 
transition  was  imposed  at  five  percent  of  the  length,  i.e.,  15.24  cm  from  the 
bow,  where  the  model  was  fitted  with  plastic  studs.  The  streamwise  velocity 
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at  the  edge  of  the  boundary  layer,  Ue,  was  determined  from  the  wave  elevation 
via  the  Bernoulli  equation 


^  i  ~  2g0g  (18) 

where  gQ  is  the  gravity  acceleration,  g  =  g(x)  is  the  wave  elevation  along  the 
hull,  and  is  the  forward  speed  of  the  ship. 

Calculations  along  other  streamlines  were  started  at  midship  with  the 
measured  velocity,  profiles  and  the  flow  was  considered  fully  turbulent.  The 
velocity  at  the  edge  of  the  boundary  layer  along  these  streamlines,  Ue>  was 
determined  from  Ug  =  u0  +  we,  where  ue  and  we  are  the  edge  values  for  velocity 
components  obtained  from  the  measurements. 

3*2.  Calculation  Results 

Figure  55  shows  the  measured  wave  profiles  on  the  hull  for  four  Froude 
numbers  of  0.266,  0.313,  0.350,  and  0.400.  The  distribution  of  the  wall  shear 
stress  presented  in  Figure  56  clearly  indicates  the  influence  of  waves.  The 
wave  system  at  Froude  numbers  of  0.266  and  0.313  are  strong  and  produce  a 
marked  oscillation  in  the  stress  distributions.  Note  that  the  usual  procedure 
of  double-model  calculations  would  result  in  a  single  unrealistic  predic¬ 
tion.  These  calculations  also  indicate  that  there  is  no  separation  at  the 
free  surface  in  the  range  of  Froude  numbers  considered.  Thus,  the  influence 
of  Froude  number  on  the  viscous  resistance  is  felt  primarily  through  the 
changes  in  the  friction  distribution  and  the  wetted-surface  area.  Figures  57 
through  59  show  the  variation  of  displacement  thickness,  momentum  thickness 
and  shape  parameter  along  the  wave  profiles  on  the  hull.  Here  again,  the 
effect  of  Froude  number  on  boundary-layer  development  is  clearly  shown  by 
stronger  oscillatory  behavior  of  the  integral  parameters  related  to  wave 
profiles  at  Froude  numbers  of  0.266  and  0.313.  Finally  Figure  60  shows  the 
variation  of  the  crossflow  angle  across  the  boundary  layer  at  three  streamwise 
positions  along  the  wave  for  Fr  =  0.313.  As  noted  earlier,  these  may  be 
regarded  as  a  measure  of  deflection  of  the  free  surface  through  the  boundary 
layer.  The  three  locations  have  been  selected  to  show  the  influence  of  the 
free-surface  curvature  which  changes  from  concave  at  ~  =  -  0.750  to  cross 
over  or  S-shape  at  =  0.380,  and  all  positive  at  ~  =  0.780.  It  is  seen 
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that  surface  deflections  are  well  correlated  with  curvature  changes.  However, 
the  deflections  are  small  since  the  Wigley  hull  is  slender  and  the  waves  are 
not  steep. 

Figure  61  shows  the  comparison  of  the  shape  parameter  along  the  first 
streamline  for  different  Froude  numbers.  Here  again  the  effect  of  Froude 
number  on  the  development  of  the  boundary  layer  is  vividly  shown  by  the  marked 
oscillation  in  the  shape-parameter  distribution  for  Fr  =  0.266  which  repre¬ 
sents  the  stronger  wave  system.  The  influence  of  the  wave  profile  on  the 
boundary  layer  lessens  toward  the  keel.  Figure  62  shows  the  variation  of  the 
shape  parameter  along  the  four  streamlines  for  Fr  =  0.266.  Note  that  the 
shape-parameter  changes  are  more  distinct  along  the  first  two  streamlines 
which  are  nearer  to  the  free  surface,  and  diminish  toward  the  keel  along 
streamlines  3  and  4.  A  similar  trend  was  noticed  in  the  results  for  the  other 
Froude  numbers,  which,  for  the  sake  of  brevity,  are  not  shown  here. 

Figure  63  shows  the  comparison  between  the  measured  and  calculated 
streamwise  momentum  thickness  for  Fr  =  0.350.  It  is  clear  that  the  calculated 
momentum  thickness  is  in  reasonable  agreement  with  the  measurements  except 
near  the  stern  where  the  boundary-layer  thickens  and  the  thin-boundary-layer 
equations  are  no  longer  accurate  enough.  The  agreement  between  calculated  and 
measured  values  is  poorer  along  the  first  streamline.  This  is  due  to  the  fact 
that  the  streamlines  begin  to  converge  toward  the  free  surface  after  mid- 
ships.  This  causes  that  boundary  layer  to  thicken  faster  near  the  free  sur¬ 
face  and  so  the  measured  values  start  to  deviate  from  the  calculated 
results.  This  trend  is  also  shown  by  the  second  streamline,  but  the  deviation 
of  the  experimental  results  from  calculations  occurs  further  downstream  since 
the  rate  of  growth  of  the  boundary  layer  along  this  streamline  is  slower  than 
for  the  first  streamline.  Finally,  along  the  fourth  streamline  near  the  keel, 
the  experimental  points  deviate  from  the  calculated  results  in  the  opposite 
direction.  This  can  be  explained  by  noting  that,  near  the  keel,  all  stream¬ 
lines  diverge  near  midships,  so  that  the  boundary  layer  is  thicker  there  and 
then  it  becomes  thinner  toward  the  stern.  The  aforementioned  is  the  general 
trend  which  was  observed  from  the  comparison  of  all  the  measured  integral 
parameters  with  their  calculated  values.  This  is  also  shown  in  Figures  64  and 
65  in  which  the  measured  shape  parameter  is  compared  with  the  calculated 
results  for  Froude  numbers  of  0.266  and  0.350.  It  is  evident  from  these 
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figures  that  the  effect  of  the  free  surface  has  been  depicted  as  much  weaker 
by  the  calculated  results  than  by  the  measurement.  This  may  be  explained  by 
noting  that  the  small-crossflow  approximation  is  in  fact  a  two-dimensional 
calculation.  Each  streamline  is  treated  separately  without  considering  the 
effect  of  neighboring  streamlines.  A  better  approach,  of  course,  would  be  a 
fully  three-dimensional  boundary-layer  calculation.  The  Wigley  hull,  however, 
is  a  slender,  or  a  'weak'  three-dimensional  body,  and  the  present  numerical 
approach  to  calculate  the  development  of  the  boundary  layer  seems  to  yield 
reasonable  results. 
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Figure  55.  Measured  wave  profiles  on  the  hull. 
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Figure  57.  Displacement  thickness  distribution  along  the  waves 
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Figure  58.  Momentum  thickness  distribution  along  the  waves. 
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Figure  59.  Shape  parameter  (H)  distribution  along  the  waves 
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Figure  60.  Variation  of  crossflow  angle  across  the  boundary 
layer  (Fr  =  0.313). 
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Figure  61.  Variation  of  shape  parameter  along  the  hull  for 
different  Froude  numbers. 
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Figure  62.  Variation  of  shape  parameter  along  different  streamlines  (Fr  =  0.266). 
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Figure  63.  Comparison  between  calculated  and  measured  momentum 
thickness  (Fr  =  0.350) . 
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Figure  64.  Comparison  between  calculated  and  measured  shape 
parameter  (Fr  =  0.350). 
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Figure  65.  Comparison  between  calculated  and  measured  shape 
parameter  (Fr  -  0.266). 


CHAPTER  V 

WAVE-RESISTANCE  CALCULATIONS 

1.  General 

In  this  chapter,  wave-resistance  coefficients  Cw  of  the  Wigley  ship  form 
will  be  calculated  with  three  different  source  distributions.  These  distribu¬ 
tions  are:  1)  Michell  thin-ship  distribution,  2)  a  first-order  source  distri¬ 
bution  on  the  hull  surface,  and  3)  a  slender-body  solution  for  the  centerplane 
distribution.  With  the  aforementioned  singularity  distributions,  we  shall 
consider  two  cases  here,  1)  wave  resistance  calculations  without  the  viscous 
effects,  and  2)  wave  resistance  with  the  effect  of  viscosity  take  into 
account.  All  the  calculations  will  be  carried  out  employing  the  Michell 
integral.  The  residuary  resistance  and  wave  resistance  obtained  from  towing 
tank  tests  will  be  used  as  references  for  comparison  of  different  distribu¬ 
tions. 


2.  Numerical  Evaluation  of  Michell (s  Integral 

for  the  Wigley  Ship  Form 

The  wave-resistance  coefficients  for  the  Wigley  hull  are  given  by  C 

12  w 
=  R/ pUqS)  where  R  is  the  wave  resistance.  The  Michell  integral  for  a  form 

with  fore-and-aft  symmetry  is  given  by 


ir/2. 


R  =  16it  p  k2  /  Q2(0  )  Sec3©  d0 


(19) 


where  Q  is  the  amplitude-distribution  function 


L/2  0 

Q  =  2  /  /  ct(x,z) 

0  -H 


exp(Kz  sec 


sin(Kx  sec0 ) 


dxdz 


(20) 


Here  a(x,z)  is  the  strength  of  the  centerplane  source  distribution,  p  the 
mass  density  of  the  fluid,  UQ  the  free-stream  velocity,  K  =  g/u2  =  (LF2)-1 
whom  g  in  the  acceleration  of  gravity,  L  is  the  ship  length,  Fr  is  the  Froude 
number,  and  0  is  the  wetted-surface  area  given  by  S  =  0.1487L^. 
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Due  to  the  oscillatory  behavior  of  the  integrands  in  (19)  and  (20),  the 
control  of  numerical  errors  in  the  calculations  is  of  primary  concern.  This 
has  been  brought  out  forcefully  by  Landweber  (1979),  who,  in  an  analysis  of 
the  computed  results  for  the  wave  resistance  of  the  Wigley  ship  form  at  a 
workshop,  showed  that  even  the  first-order  calculations  produced  discrepancies 
of  the  order  of  15  percent.  In  the  following,  a  procedure  in  which  the  oscil¬ 
latory  part  of  the  Integrand  is  integrated  exactly  has  been  developed. 

X-AX/Z  X+  4X/£ 

Consider  a  small  element  of  the  A^/Z 

ship  centerplane  with  dimensions 
Ax  and  Az  in  the  x  and  z 

directions,  respectively.  Let  us  *--=£■  AZ/t 

Z  -  AZ/Z 

assume  that  the  strength  of  the  source  distribution  over  this  element  can  be 
approximated  by  cr_  which  is  constant  over  the  area  AxAz.  Thus  the  integral 
(20)  for  this  element  can  be  written  as 

_  x  +  Ax/2  z  +  Az/2 

Qj  3  =  2a4 ,  /  Sin  (Kx  Sec9  )  dx  /  exp  (Kz  Sec^Q  )  dz 


ij 


ij 


x  -  Ax/2 


z  -  Az/2 


(21) 


which  is  exactly  integrable  and  yields 

8a .  .  .  o 

Q.  .  =  ■x  Sin  (K  4s  Sec  9 )  Sinh  (K  Sec  9 )  Sin  (Kx  .  Sec9  ) 

1J  K  Sec^Q 

•  exp  (Kz.  Sec29  )  (22) 

J 

After  summation  over  the  entire  centerplane  and  substitution  in  (19),  we 
obtain 

R  =  lOS^rpU2  j  E(0)  co2  e  d9  (23) 

YT  o 

where 

l^EO)  =  sin2(K  xr  Sec9  )  sinh2(K  Sec29  ) 

M1  M2 

(  l  \  a..  Sin  ( Kx .  Sec9  )  exp  (Kz.  Sec29  )) 2  (  24) 

i=l  j=l  1J  1  J 
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and  M-|_  and  M2  are  numbers  of  step  sizes  in  x  and  z  directions,  respectively. 


Integration  by  parts  in  (23)  yields 
ir/2 

E(0 )  Cos^e  d0  =  -  j 

o 


o  IT  /2  . 

/  E(0 )  Cos30  d0  =  -  jr  /  ( +  3  sin  0)  ^  d0 


dE 


(25) 


The  range  of  integration  is  then  subdivided  into  N  equal  subintervals 
at  0  =  0^,  0£,  =  \  and  the  quadrature  formula 


0 

p-q 


1  ,sln39i*i 


0 


i+l,sin30  .  „  .  A,  dE  ,A  -  A  , 

' — 3 —  +  3  sm0  )  -gg-  d0  =2  ( - 3 

sin30  . 

. 3 -  1  +  3  sin91)(Ei+1  -  E.) 


+  3  sin0i+1  + 


(26) 


by  which  the  rapidly  varying  factor  E  is  treated  exactly,  is  assumed. 
Here  Ei  =  E(0^).  Substituting  (26)  and  (25)  into  (23),  we  obtain 

,0{J  tt2  N  sin30  sin  30 

R  '  <  En  En+1  - f1  -  En  +  ^En+lsin9n+l  ' 

iv  n-i 

sin30  sin30 

+  En+i  — +  3En+lsin0n  "  En  - “  3E~sln0  —  5  (27) 


Since 


N  sin 30  sin30 

EE  n+1  -  n 


sin30 


n=l 

and  similarly 


!n  — 1 - -  \  — 


=  E 


N+l 


is  '  X 

'N+X  3  -  Eq  sin30-^  =  0 


N 

3(Enil  sinen+1  -  Ensinfin)  =  3(ENtlsi„9Ntl  -  E-pinep  =  0 

We  see  that  the  sum  of  the  first  four  terms  in  (27)  is  zero.  The  remaining 
terms  then  yield 


0<-,  t2  N-l  sin30 

P  _  -256iL  ^4  r  r  w  , _ _ 

Cw  S  Fr  ^  En+1  (  3 

n=l 


n 


+  3sin0  )  -  E 
n'  n 


sin  30 


n+l  .  _ _  .  8 


~3 -  +  3sin0n+l)  -  3  V 


(28) 


where  E(0^)  is  given  by  (24). 
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The  accuracy  of  the  integration  with  respect  to  x,  z  and  9  was  investi¬ 
gated  by  halving  the  step  size  for  each  variable  while  keeping  the  other  two 
fixed.  By  this  procedure,  step  sizes  Ax  =  -^j,  Az  =  and  A0  =  were 
found  to  give  sufficiently  small  errors. 


3.  Wave  Resistance  without  the  Effects  of  Boundary  Layer 

In  this  part,  the  wave  resistance  of  the  Wigley  form  was  calculated  for 
several  different  source  distributions,  on  the  centerplane  or  on  the  hull, 
without  considering  the  effects  of  the  boundary  layer  on  the  wavemaking  resis¬ 
tance.  As  was  mentioned  earlier,  all  the  wave  resistance  calculations  are 
carried  out  only  to  the  first  order. 


3.1.  Thin-Ship  Approximation 

The  Michell  first-order  approximation  for  the  source  distribution  for  the 
Wigley  form  is  given  by 


a (x,z) 


U  f  2 

o  x  _  -Bx  ,  z 

2tt  Dr  u  "  2 

n 


(29) 


where  y  =  f(x,z)  is  the  equation  of  the  hull  surface  and  fv  =  With  this 
distribution,  (20)  is  exactly  integrable  and  the  resulting  equation  is  in  much 
simpler  form  than  (28);  see  Landweber  (1979).  Here  (29)  was  treated  as  a 
general  distribution  and  substituted  into  (24)  and  (28)  for  calculation  of  the 
wave-resistance  coefficients.  The  results  obtained  were  exactly  the  same  as 
those  obtained  by  Landweber  and  Celik  (1980).  Figure  66  shows  the  comparison 
between  this  calculation  and  the  wave-resistance  coefficient,  C  =  (\  -  C 

7  WV  T  V> 

obtained  from  towing-tank  measurements.  Here  Ct  is  the  total  resistance 
coefficient,  Cv  the  viscous  resistance  coefficient  obtained  by  wake  survey, 
and  is  the  measured  wave-resistance  coefficient.  This  figure  also  shows 
the  comparison  between  the  calculation  and  residuary  resistance  C„  =  C+  - 
Citto  where  Cj^q  is  the  International  Towing  Tank  Committee  (ITTC)  1957 
friction  line  given  by  CITTC  =  0.075  (log10  Re  -  2)“2. 


3.2.  Hull-Surface  Distribution 

An  alternative  first-order  source  distribution  on  the  hull  surface, 
derived  as  the  first  approximation  to  the  solution  of  the  Fredholm  integral 
equation  of  the  second  kind  at  zero  Froude  number,  is  given  by 
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4 -no  (£  ,?  )  = 


(30) 


where  n  =  ±  f(£,?)  is  the  equation  of  the  Wigley  hull  surface,  D  =  (1 

2  2  1/2 

+  f^  +  f  )  7  ,  and  £  and  c  are  coordinates  of  the  points  on  the  hull.  A.  fur¬ 
ther  improvement  is  given  by  the  added-mass  correction 


o (x,z) 


:i+kn  )U  f 
1  ox 

4rr  D 


(3D 


Here  the  longitudinal  added-mass  coefficient  k-^  can  be  estimated  from  its 
known  value  for  an  'equivalent  ellipsoid' ,  one  having  the  same  volume,  mid¬ 
section  area  A^,  and  free-surface  area  Afg  as  the  original  ship  hull.  For 
the  Wigley  ship  form,  V  =  BLH,  A^  =  |  BH,  and  Afs  =  |  BL,  where  L,  H  and  B 
are  length,  draft,  and  beam  of  the  model,  respectively.  The  equivalent  ellip¬ 
soid  with  these  sectional  areas  and  volume  has  an  axis  ratio  of  8:  vf-  it  •  1 

15 

(Miloh  and  Landweber,  1980),  and  the  added-mass  coefficient  for  this  ellipsoid 
is  k^  =  0.0240. 

Since  the  source  distribution  in  (31)  is  a  surface  distribution  rather 
than  a  centerplane  distribution,  equations  (19)  and  (20)  need  to  be  rede¬ 
rived.  The  wave  potential  for  a  surface  distribution  of  sources 
with  n  =  f(£,?)  is  given  by:  (Wehausen  and  Laitore,  I960) 

*/2 

<l>(x,y,z)  =  -  4K  ff  a(£,s)  {cos  [K(y-n)sin0Sec  0]  +  cos[K(y+n) 
s  o 

(32) 

sin0Sec20]}  Sin[K(x-£  )Sec0  ]exp[K(z+c  )Sec20  ]Sec20d0dS  +  <j>2(x,y,z) 

By  applying  the  Lagally  theorem,  one  can  calculate  the  resistance  due  to  this 
source  distribution, 

R  =  -  4ttp  /  o  (x,z)  dS  (33) 

s 

It  should  be  noted  that  the  <j>2(x,y,z)  in  (32),  which  is  a  (K,0  )  double  inte¬ 
gral,  does  not  contribute  to  the  Lagally-computed  force  (Lunde,  1951),  and 
therefore  has  not  been  written  in  detail  in  (32)  and  will  be  omitted  in  the 
following  calculations. 
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From  (32),  the  induced  velocity  can  be  written  as 


9<j> 

5x  " 


2  L/2  0  tt/2 

-  8K  /  /  /  c(5,c)  cos  (KysinGSec  0)  cos  (Kr\sin0Sec  0  ) 

-L/2  -H  0 

cos  [K(x-£  )Sec0  ]  exp  [K(z+£  )Sec20  ]Sec^0D(£  )dg  d?  d0 


(34) 


Here,  the  element  of  the  hull  surface  has  been  projected  onto  the  centerplane, 
i.e.,  dS  =  Dd£d?  where  D(£ ,? )  =  (1+  f2  +  f2)1/2.  Substituting  (34)  into  (33) 
and  considering  the  source  distribution  for  both  sides  of  the  hull,  we  obtain 

7T* 

2 

R  =  64ffpK2  /  Q2Sec^0d0  (35) 

o 

where  Q  is  the  amplitude-distribution  function 


L 

2 


Q  «/ 


0 

/  a(x,z)Sin(KxSec0 )cos(KySin0Sec20 )exp(Kzsec20 )D(x,z)dxdz 

-H 


(36) 


Here  y  =  f(x,z)  is  the  equation  of  the  Wigley  hull. 

Wave-resistance  coefficients  Cw  for  the  surface  distribution  (31)  were 
calculated  by  using  relations  (35)  and  (36).  Figure  67  shows  the  results  of 
this  calculation  and  its  comparison  with  the  Michell  thin-ship 
approximation.  The  comparison  shows  that  the  results  with  the  hull-surface 
distribution  (30)  are  in  poorer  agreement  with  the  residuary  Cr  and  measured 
wave  resistance  Cw  than  those  from  Michell  formula  (29). 


3.3*  Slender-Body  Solution  for  the  Centerplane  Distribution 

Landweber  and  Ju  (1983),  obtained  a  slender-body  solution  in  terms  of  a 
centerplane  source  distribution  for  a  'modified'  Wigley  ship  form.  Their 
modification  is  that  some  transverse  sections  of  the  model  were  slightly 
deformed.  Ship  sections  for--j<x<-j  are  gradually  varied  from  a  parabolic 
section  in  such  a  way  that  the  midship  section  is  transformed  into  an  ogive 
with  the  same  section  area.  This  centerplane  source  distribution  was  used  in 
(19)  for  calculation  of  the  wave-resistance  coefficients.  Figure  68  shows  the 
result  of  this  calculation  and  its  comparison  with  the  Michell  distribution 


106 


FROUDE  NUMBER 

Figure  66  Comparison  between  Michel!  integral  and  residuary 
resistance. 
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FROUDE  NUMBER 

Figure  67.  Comparison  between  hull-surface  distribution  and 
residuary  resistance. 
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Figure  68.  Comparison  between  slender-body  solution  and 
residuary  resistance. 
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and  the  residuary  resistance.  The  agreement  between  residuary  resistance  and 
this  calculation  is  very  reasonable  up  to  about  Fr  =  0.30,  and  then  deviates 
sharply.  This  may  be  expected  since  the  slender-body  distribution  as  well  as 
the  thin-ship  distribution  was  obtained  for  the  zero-Froude-number  case  and  it 
cannot  be  expected  to  be  accurate  at  high  Froude  numbers. 


4»  Wave  Resistance  With  the  Effect  of  Boundary  Layer 

The  relationship  between  the  flow  exterior  to  a  boundary  layer  and  wake, 
and  the  concept  of  displacement  thickness  and  the  source  distribution  which 
generates  the  outer  irrotational  flow  were  examined  in  detail  for  the  case  of 
two-dimensional  and  axisymmetric  bodies  by  Preston  (1945),  Lighthill  (1958), 
and  Landweber  (1978).  In  continuation  of  his  work,  Landweber  (1981),  using  an 
extension  of  the  analysis  used  in  the  two-dimensional  case,  derived  a  second- 
order  approximation  for  a  centerplane  source  distribution  for  a  ship  form. 
The  first-order  expression  for  the  centerplane  distribution  of  sources  equiva¬ 
lent  to  the  boundary  layer  of  a  ship  form  is  given  by: 
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a ( x,  z )  =  o  (x,z)  + 

x>' 


2tt  £ 
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where  a^(x,z)  is  the  potential-flow  distribution  and  the  rest  of  the  terms  are 

due  to  viscous  corrections.  {£ .  .}  denotes  the  matrix  of  direction  cosines 

1J 

relating  the  (x,y,z)  coordinates  of  the  ship  form  to  a  right-handed,  curvilin¬ 
ear,  surface-related,  nonorthogonal  coordinate  system  (s,t,n),  where  n  denotes 
distance  in  the  direction  of  the  outward  normal  to  the  hull  surface  S,  and  n  = 
0  on  S.  Here  s  =  constant  in  the  zy  plane  and  t  -  constant  in  the  xy  plane. 
The  unit  vectors  for  the  (s,t,n)  coordinate  system  are  related  to  the  ship 
hull  equation  y  =  f(x,z),  by 
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where  (i,j,k)  are  unit  vectors  in  the  cartesian  coordinate  system  (x,y,z). 

Note  that,  for  the  (s,t,n)  coordinate  system,  e^  »~e  =  e^  • "e  =  0 

s  n  t  n 

but  eg*et  *  0.  The  latter  product,  which  is  a  measure  of  the  nonorthogonality 
of  the  coordinate  system,  obtained  from  (38),  is 


e  = 
s 


f  f 
x  z 


(1  +  f?)  (1  +  f2) ]1/2 

X  Zi 


(39) 


and  has  its  maximum  e^  eg  =  0  at  midship  and  stem,  and  its  minimum 


e  •  e  = 
t  s 


-  0.041  at  a  distance  about  L/3  from  the  midship,  and  H/2  from  the  keel, 
and  K2  in  (37)  are  the  curvatures  of  surfaces  s  =  const,  and  t  =  const., 
respectively,  and  since  these  surfaces  are  the  zy-  and  xy -planes,  then 


K, 


=  0. 


Definitions  of  6  1  and  6 g  are  very  similar  to  that  of  the  displacement 
thickness, 

6  ug-  u  6  w  -  w 

61  =  /  “u -  dn,  «2  =  /  -  dn  (40) 

O  00  o  00 

where  (u,w,v)  are  the  velocity  components  in  the  (s,t,n)  directions,  and  ue> 
we  denote  values  of  u,w  at  n  =  6 . 

As  was  mentioned  in  Chapter  II,  the  boundary-layer  measurements  were 
originally  taken  in  constant  waterline  planes.  In  order  to  calculate  the 
various  terms  in  (37),  however,  we  need  to  have  velocity  components  in  the 
(s,t,n)  coordinate  system.  To  solve  this  problem  a  software  was  developed  in 
which  first  a  linear  interpolation  between  the  measured  velocity  components  on 
each  profile  was  employed,  and  then  these  profiles  were  interpolated  in  the  x 
and  z  directions  by  least-square  polynomials.  The  whole  domain,  was  then 
divided  into  small  cubical  elements  and  the  three  velocity  components  were 
calculated  from  interpolating  polynomials  at  their  eight  nodal  points.  The 
coordinates  of  points  along  normals  to  the  hull  were  calculated  from  the 
geometry  of  the  Wigley  hull.  The  velocity  components  at  each  point  on  these 
normals  were  calculated  by  linear  interpolation  among  the  eight  nodal  points 


ill 


of  the  cubical  element  surrounding  the  desired  point.  The  obtained  velocity 
components  on  the  normals  to  the  hull  (u,v,w)  were  along  cartesian  coordinates 
(x,y,z),  respectively.  The  following  relation  was  employed  to  transfer  them 
to  (us,  u-j.)  along  (s,n,t)  coordinate  lines: 

t  =  ui  +  vj  +  wk  =  ugeg  +  un?n  +  utet  (41) 
where  (e  ,  e  ,  e.)  are  unit  vectors  given  by  (38). 

o  II  U 

The  complete  sets  of  (ug,  ut,  un)  for  Froude  numbers  0.266,  0.313,  0.350, 
and  0.400  are  stored  on  magnetic  tapes  and  are  available  from  IIHR  tape 
library.  See  Appendix  A. 

Using  these  velocity  profiles,  the  expressions  (40)  for  ^  and  6  were 
calculated,  (see  Tables  5. a  to  5.d)  and  then,  employing  <5^  and  6^  and  their 
derivation,  the  viscous-correction  terms  of  equation  (37)  were  calculated  for 
all  four  Froude  numbers  and  by  interpolation,  at  several  other  Froude  num¬ 
bers.  The  result  was  a  set  of  centerplane  source  distributions  which  extended 
from  midship  to  the  stern.  The  wave-resistance  calculations  for  the  three 
aforementioned  cases  were  repeated,  using  expression  (37)  as  the  source  dis¬ 
tribution  term  in  the  appropriate  equations.  Figure  69  shows  the  comparison 
with  the  Miehell  thin-ship  source  distribution  with  and  without  the  viscous- 
correction  term. 

4.1.  Miehell  Source  Distribution  with  Viscous  Effects 

With  equation  (37)  as  the  centerplane  source  distribution  in  the  Miehell 
integral,  the  wave-resistance  coefficients  for  the  four  Froude  numbers  0.266, 
0.313,  0.350,  and  0.400  were  calculated.  Figure  70  shows  the  comparison 
between  Cw,  Cr  and  Miehell  centerplane  distribution  with  and  without  vis¬ 
cosity  effects.  It  is  obvious  that  including  the  viscosity  effects  in  wave- 
resistance  calculations  modifies  the  theoretical  resistance  in  the  right 
direction,  e.g.,  at  Fr  =  0.313  this  lowers  the  curve,  whereas  at  Fr  =  0.350 
this  elevates  the  curve  so  that  it  is  in  much  better  agreement  with  the  exper¬ 
imental  results. 
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Table  5.  Displacement  thickness  of  Eq.  40. 
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Table  5.  Continued 
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Figure  M 


Comparison  between  cr(x,z)  and  <s 
Michel!  canterplane  distribution 


(x,2),(Eqv37),for 
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4.2.  Slender-Body  Solution  with  Viscous  Effects 


With  the  slender-body  solution  for  a^(x,z)  in  equation  (3V),  the  wave- 
resistance  coefficients  for  the  modified  Wigley  form  were  calculated.  Figure 
71  shows  the  comparison  between  CW)  residuary  resistance,  and  the  slender- 
body  solution  with  and  without  viscosity  effects.  Here  again  it  is  apparent 
that  by  talcing  effect  of  boundary  layer  into  account,  the  theoretical  results 
improve . 


4.3.  Hull-Surface  Distribution  with  Viscous  Effects 

In  order  to  take  into  consideration  the  effect  of  the  boundary  layer  on 
wave  resistance  for  hull-surface  distribution,  we  had  to  modify  the  Michell 
integral.  This  is  because  the  viscosity  effect  is  in  the  form  of  a  center- 
plane  distribution  of  singularities  while  the  potential  flow  is  in  the  form  of 
a  hull-surface  distribution.  Let  subscript  H  and  CP  refer  to  hull  and  center- 
plane,  respectively.  The  potentials  due  to  hull  and  centerplane  distribu¬ 
tions,  neglecting  the  free-wave  terms  are 

L/2  0  n/2 

(f)„(x,y,z)  =  -  8K  f  f  f  exp[K(z-c  )Sec  9  ]Sin[K(x-?  )Sec9  ] 

-L/2  -H  0  (42) 

Sec20aH(£  ,?  )Cos(KySin0Sec20  )Cos(KnSin0Sec20  )D(£  ,£  )d£  d?  d0 


L/2  0  tt/2 

^rp(x, y,z)  =  -  8K  /  /  /  exp[K(z-^  )Sec  0  ]Sin[K(x-£  )Sec0  ] 

-L/2  -H  0 

Sec20agp(£  ,?  )Cos(KySin0Sec20  )dC  d£  d0 


(43) 


where  y  =  f(x,z)  and  n  =  f (£ ,C )  are  equations  of  the  hull.  The  Lagally  theor¬ 
em  (33),  gives  two  force  components,  and  R2,  due  to  the  surface  and  center- 
plane  distributions,  respectively, 


R,  =  -  8irp 
R2  =  -  4irp 


L/2 

/ 

-L/2 

L/2 

/ 

-L/2 


0  9<f>H  3<f>rp 

/  [  +  (~3Y”)H'aH^x,z^D^x,z^dxclz 

■H 

0  9<f>H  a<f,rp 

l  ^jTx^cp  +  (‘1T~)cp]  acp  (x,z)  dx  dz 

•H 


(44) 

(45) 
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Differentiating  (42)  and  (43)  with  respect  to  x  and  substituting  into  (44)  and 
(45)  gives  the  total  resistance  R  =  R1  +  r2, 

K 

o  2 

R  =  32  tt pIT  /  (2P+Po)(P+Po)  Sec30  d9  (46) 

where 


L/2  0 

P(9 )  =  /  /  a  (x,z) 

-L/2  -H  H 


exp(KzSec20  )Sin(KxSec9  )Cos(KySin9Sec29  ) 


D(x,z)  dxdz 

L/2  0 

P  (9)  =  /  f  a  (x, z)  exp(KzSec29 )Sin(KxSec9 )  dxdz 
-L/2  -H 


(47) 

(48) 


The  wave-resistance  coefficients  were  calculated  by  using  equation  (46). 
Figure  (72)  shows  the  comparison  among  the  residuary  resistance  coefficient, 
^wv>  and  calculated  from  the  hull-surface  distribution,  with  and  without 
boundary  layer  effects.  Here  again  the  theoretical  results  are  in  much  better 

agreement  with  residuary  resistance  and  Gm  when  the  boundary-layer  effects 
are  taken  into  account. 
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FROUDE  NUMBER 

Figure  70.  Comparison  between  Michell's  thin  ship  with  and 
without  viscous  affects. 
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FROUDE  NUMBER  (FR) 

Figure  7^  Comparison  between  slender-body  solutions  with 
and  without  viscous  effects. 
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FROUDE  NUMBER 

Figure  72.  Comparison  between  hull-surface  distributions  wi 
and  without  viscous  effects. 
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5.  Discussion  of  the  Results 


The  wave-resistance  calculations  for  the  Wigley  ship  form  without  consid¬ 
ering  the  effect  of  viscosity  showed  that  the  Michell  thin-ship  theory  yields 
better  results  than  the  first-order  hull  distribution.  A  similar  result  has 
been  pointed  out  by  Breslin  (1963),  and  Kang  (1978).  It  may  be  reasoned  that 
the  second-order  terms  due  to  the  nonlinearity  of  the  hull  boundary  condition 
and  those  due  to  the  nonlinearity  of  the  free-surface  condition  approximately 
compensate  each  other  in  the  Michell  integral,  which  gives  the  wave  resistance 
from  the  first-order,  fully-linearized  theory  for  a  thin  body.  In  the  hull- 
surface  distribution,  however,  only  the  hull  boundary  condition  has  been 
improved.  This  indicates  that  all  the  second— order  terms  in  the  second— order 
thin-ship  theory  are  quite  important  in  wave-resistance  calculations  and  one 
should  consider  all  of  them  to  get  any  better  results  than  the  Michell  thin- 
ship  distribution. 

Table  6  shows  a  comparison  of  various  results  obtained  in  this  chapter. 
Including  the  effect  of  viscosity  in  a  wave-resistance  calculation,  even  if  it 
is  limited  to  first-order  terms  only ,  results  in  a  great  improvement  in  the 
Michell  and  hull-surface-distribution  calculations.  It  reduces  the  percentage 
of  the  difference  between  theory  and  experiment  by  more  than  50$.  In  the  case 
of  the  modified-slender-body  solution,  however,  including  the  viscous  effects 
does  not  appear  to  change  the  results  by  much.  This  may  be  due  to  the  fact 
that  viscous  corrections  are  for  a  Wigley  hull  with  parabolic  section  and  not 
the  modified  form.  Even  in  this  case,  however,  introducing  the  viscous 
effects  reduces  the  peaks  of  the  Michell  integral,  resulting  to  a  much  smooth¬ 
er  curve. 

It  should  be  noted  that,  in  this  section,  contrary  to  Kang's  (1978) 
analysis,  we  did  not  include  the  effect  of  wake  sources  on  the  wavemaking 
resistance  of  the  Wigley  hull.  This  is  due  to  the  fact  that  these  sources, 
which  represent  the  displacement  thickness  of  the  wake  and  are  known  as  Betz 
(1925)  sources,  contribute  mainly  to  the  viscous  drag  of  the  ship  model  and 
their  effect  on  wave  resistance  is  negligible.  It  has  been  shown  experimen¬ 
tally  by  Swain  &  Landweber  (1976)  that  the  wake  of  a  ship  model  would  contri¬ 
bute  only  about  one  percent  of  the  total  wave  resistance. 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  most  useful  purpose  served  by  the  present  study  is  the  collection  of 
a  complete  set  of  data  on  boundary-layer,  pressure-distribution,  and  resis¬ 
tance  measurements  for  the  Wigley  parabolic  hull.  These  data,  to  the  best  of 
our  knowledge,  are  the  first  complete  set  of  boundary— layer  measurements  which 
have  been  obtained  for  a  ship  model  at  zero  trim  and  sinkage.  This  is  very 
significant,  because  most  computation  methods  are  made  for  the  fixed-model 
condition.  For  example,  if  one  cannot  predict  wave  resistance  for  the  fixed- 
model  condition,  a  good  prediction  for  the  more  realistic  f ree-to-trim-and- 
sink  condition  is  unlikely,  since  one  has  to  use  the  fixed-model  condition  as 
the  initial  condition  of  an  iterative  procedure.  In  short,  the  data  obtained 
in  the  present  study  can  serve  as  the  basis  for  further  theoretical  studies  of 
the  Wigley  hull.  The  major  conclusions  of  this  study  are  as  follows: 

1)  Experimental  results  for  the  total-resistance  coefficients  of  the 
Wigley  parabolic  hull  for  free  and  fixed  conditions  are  shown  in  Tables  1  and 
2,  respectively.  The  comparison  between  these  two  tables  shows  that  the  hump 
and  hollows  for  both  cases  occur  at  the  same  Froude  numbers,  but  the  values 
with  the  free  condition  are  greater  than  those  with  the  fixed  condition  with 
the  increment  increasing  to  about  ten  percent  at  the  largest  Froude  number  of 
0.400. 

2)  Experimental  results  for  the  pressure-distribution  coefficients  Cp  of 
the  Wigley  hull,  restrained  in  both  trim  and  sinkage,  are  shown  in  Table  3. 
An  important  feature  of  these  measurements  is  the  fact  that  the  pressure 
distribution  on  the  hull  closely  follows  the  trend  of  the  wave  profile  at  the 
free  surface,  i.e.,  shows  the  same  hollows  and  humps.  This  shows  that  the 

Froude  number  plays  an  important  role  on  the  boundary  layer  and  its  develop¬ 
ment  on  ship  hulls. 

3)  Three-dimensional  boundary-layer  measurements  were  carried  out  for 
four  Froude  numbers  of  0.266,  0.313,  0.350  and  0.400.  The  three  components  of 
the  velocity  vector  were  obtained  at  ten  transverse  sections  and  five  water- 
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lines  at  each  section,  from  midship  to  the  stem.  These  measurements  were 
extended  into  the  near  wake  of  the  model  up  to  a  distance  of  0.1  model  lengths 
downstream  from  the  stern.  The  wake  measurements  were  carried  out  at  four 
sections  and  five  waterlines  at  each  section.  The  results  of  measurement  are 
given  in  Appendix  B. 

The  boundary -layer  and  wake  measurements  show  no  flow  reversal  or  separ¬ 
ation  for  the  Wigley  hull.  This  indicates  that  this  model  is  a  very  good 
research  vehicle  for  comparison  of  the  analytical  results  or  numerical  pro¬ 
cedures  with  experimental  data. 

4)  The  boundary-layer  measurements  showed  that  the  flow  over  the  Wigley 
hull  may  be  regarded  as  one  of  small  crossflow,  since  the  measured  crossflow 
angles  are  less  than  ten  degrees. 

5)  The  original  measurements,  which  were  taken  in  transverse  sections 
and  along  constant  waterlines,  were  transferred  by  interpolation  to  data  along 
normals  to  the  hull. 

6)  Boundary-layer  calculations  on  the  Wigley  hull  were  performed  using 
the  small  crossflow  approximation  of  the  boundary-layer  equations.  These 
calculations  were  carried  out  along  the  streamlines  on  the  hull  obtained  using 
Guilloton's  method. 

The  comparison  of  the  calculated  and  measured  results  show  good  agreement 
over  most  of  the  hull  except  near  the  stem  where  the  boundary  layer  is  rather 
thick. 

7)  Wave-resistance  coefficients  of  the  Wigley  hull  were  calculated  by 
using  the  Michell  integral.  These  calculations  were  carried  out  by  employing 
the  Michell  thin-ship  centerplane  distribution,  a  first-order  hull-surface 
distribution,  and  a  slender-body  source  distribution.  Comparison  of  the 
results  showed  that  the  Michell  thin-ship  distribution  gives  a  better  estimate 
of  the  wave-resistance  coefficients  than  the  hull-surface  distribution,  but 
the  gold  medal  should  go  to  the  slender-body  solution  which  gives  the  best 
comparison  with  the  measured  wave  resistance  up  to  Fr  =  0.3.  The  agreement 
between  experimental  results  and  all  the  aforementioned  calculations  becomes 
poorer  with  increasing  Froude  number.  This  is  due  to  the  fact  that  the  calcu¬ 
lations  are  based  on  linearized  theory  and  zero-Froude-number  distributions 
and  cannot  be  expected  to  yield  satisfactory  results  at  high  Froude  numbers. 
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8)  The  equivalent  irrotational  flow  exterior  to  the  boundary  layer  was 
obtained  by  generating  the  potential  flow  exterior  to  the  boundary  layer  by 
means  of  a  centerplane  source  distribution  which  satisfies  the  boundary  condi¬ 
tion  on  the  edge  of  the  boundary  layer. 

9)  Using  the  equivalent  irrotational  flow,  the  wave-resistance  coef¬ 
ficients  of  the  Wigley  ship  form  were  calculated  in  combination  with  the 
distributions  mentioned  in  (7)  to  study  the  effect  of  viscosity  on  the  wave¬ 
making  resistance  of  the  ship  model.  These  calculations  showed  that  the 
agreement  with  the  "measured"  resistances  was  considerably  improved  by  modify¬ 
ing  the  source  strength  for  the  effects  of  viscosity. 

10)  The  comparison  between  different  methods  of  calculation  showed  that 
introducing  the  corrections  for  viscous  effects  improved  the  results  with  the 
Michell  thin-ship  and  the  hull-surface  distribution  more  than  with  the  slen- 
der-body  solution  although  the  latter  needed  only  a  small  correction. 

All  the  wave-resistance  calculations  in  the  present  study  were  limited  to 
first-order  terms.  For  future  work,  it  is  recommended  that  second  and  higher- 
order  terms  be  taken  into  account.  The  slender-body  solution,  in  particular, 
may  be  the  best  answer. 

Although  the  effects  of  viscosity,  in  the  present  study,  were  confined 
only  to  the  boundary  layer  over  half  of  the  ship  model,  the  results  showed 
marked  improvement.  A  better  agreement  between  theory  and  experiments  may  be 
obtained  if  one  considers  the  boundary  layer  and  wake  of  the  whole  model,  and 
since  boundary-layer  measurements  in  a  towing  tank  are  formidable  task,  the 
next  best  alternative  would  be  to  calculate  the  ship  boundary  layers  numer¬ 
ically,  for  nonzero  Froude  numbers.  Development  of  a  reliable  method  for 
nonzero  Froude-number  calculations  of  a  ship's  boundary  layer  should  have  high 
priority. 
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